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Abstract 


Objectives  of  this  study  were  to  measure  the  extent  of 
changes  in  organic  carbon  in  cultivated  as  compared  to 
native  soils  of  all  climatic  zones  and  several  cropping 
systems  throughout  Alberta;  and  to  examine  the 
interrelationships  among  observed  changes  in  soil  organic 
carbon  content,  soil  properties  and  regional  or  management 
variables . 

Profile  samples  were  collected  from  cultivated  and 
adjacent  native  sites  on  72  farms  located  in  the  Brown,  Dark 
Brown,  Black,  Dark  Gray  and  Gray  soil  zones  of  southern  and 
central  Alberta  and  from  the  Peace  River  area  together  with 
data  on  management  and  site  history. 

Concentration  losses  of  organic  carbon  in  A  horizons 
ranged  from  13%  to  75%  while  mass  losses  ranged  between  15% 
and  31%.  For  total  profiles  the  mass  loss  was  between  21% 
and  39%.  Variable  gains  in  mass  of  organic  carbon  in  A 
horizons  occurred  on  13%  to  43%  of  the  farms  depending  on 
the  sampling  zone. 

There  was  no  consistent  decline  of  organic  carbon  with 
time.  The  change  in  organic  carbon  was  related  to  management 
practices  and  original  nature  of  the  soil.  Losses  in  organic 
carbon  tended  to  be  lower  under  continuous  cropping  than 
under  summer  fallow . 

Changes  in  management  from  fallowing  to  more  continuous 
cropping  seem  to  have  reversed  the  direction  of  changes  in 
organic  carbon.  Existing  models  on  organic  matter  dynamics 


IV 


, 


will  have  to  be  revised  to  accommodate  changes  in 
management . 


- 


Acknowledgements 

I  wish  to  express  my  grateful  and  sincere  appreciation  to 
the  following  people  and  institutions: 

The  Alberta  Heritage  Trust  Fund  for  financial  support  of  the 
project . 

The  German  Academic  Exchange  Association  (DAAD)  for  academic 
support . 

My  supervisor  Dr.  McGill  for  his  time,  constructive 
criticism  and  guidance. 

My  committee  members  Dr.  Pettapiece,  Dr.  Robertson  and  Dr. 
Taerum  for  their  encouragement  and  advice. 

My  field  assistants  Josie  van  Lendt  and  Lucy  van  Esbroeck 
for  dedicated  help. 

All  the  persons  who  were  involved  in  tedious  sample 
preparation  and  lab  analyses. 

The  Alberta  Agriculture  Plant  Industry  Personnel,  Regional 
Soils  Specialists  and  District  Agriculturalists  for 
important  assistance  in  the  project. 

All  the  farmers  for  their  time,  cooperation  and 
encouragement  which  was  invaluable. 

My  husband  Tim  for  his,  patience,  understanding  and  moral 
support . 


vi 


Table  of  Contents 

Chapter  Page 

1 .  Introduction  . 1 

2.  Literature  Review  . 4 

2.1  Introduction  . 4 

2.2  Organic  Matter  Dynamics  . 5 

2.3  Effects  of  Cultivation  on  Organic  Matter  . 12 

2.4  Conclusions  . 31 

3.  Materials  and  Methods  . 33 

3  .  1  Field  Work  . 33 

3.1.1  Site  Selection  . 33 

3.1.2  Soil  Sampling  . 35 

3.1.3  Field  History  . 37 

3.2  Lab  Procedures  . 38 

3.2.1  Sample  Preparation  . 38 

3.2.2  Apparent  Bulk  Density  . 38 

3.2.3  Soil  pH  . 39 

3.2.4  Inorganic  Carbon  . 39 

3.2.5  Total  Carbon  . 41 

3.2.6  Organic  Carbon  . 42 

3.2.7  Moisture  Retention  . 42 

3.3  Data  Management  . 43 

4.  Changes  in  Soil  Properties  . 44 

4.1  Introduction  . 44 

4.2  Methods  of  Statistical  Analysis  . 44 

4.3  Results  . 46 


4.3.1  Calculation  of  Changes  in  Organic  Carbon  ..46 


VI  1 


4.3.2  Concentration  of  Organic  Carbon  in 

Cultivated  and  Native  A  Horizons  . 49 

4.3.3  Mass  of  Organic  Carbon  in  Cultivated  and 

Native  A  Horizons  and  Profiles  . 51 

4.3.4  Changes  in  Concentration  and  Mass  of 

Organic  Carbon  in  A  Horizons  . 53 

4.3.5  Average  Thickness  of  Cultivated  and 

Native  A  Horizons  . 53 

4.3.6  Average  Apparent  Bulk  Density  of 

Cultivated  and  Native  A  Horizons  . 56 

4.3.7  Average  pH  of  Cultivated  and  Native  A 

Horizons  . 56 

4.3.8  Distribution  of  Organic  Carbon,  Apparent 

Bulk  Density  and  pH  with  Depth  in 
Cultivated  and  Native  Profiles  . 56 

4.3.9  Moisture  Retention  in  Cultivated  and 

Native  A  Horizons  . 65 

4.4  Discussion  . ....67 

4.5  Conclusions  . 78 

5.  Regional  and  Management  Variables  Affecting  Changes 

in  Soil  Organic  Carbon  . 81 

5.1  Introduction  . 81 

5.2  Methods  of  Statistical  Analysis  . 81 

5.3  Results  . 82 

5.3.1  Characteristics  of  Management  Practices 

throughout  Alberta  . 82 

5.3.2  Relationship  between  Crop  Rotation  and 

Change  in  Organic  Carbon  . 86 

5.3.3  Simple  Linear  Regression  Models  . 86 

5.3.4  Multiple  Linear  Regression  Models  . 101 

5.4  Discussion  . 112 

5.5  Conclusions  . 117 

6.  Overall  Summary  and  Conclusions  . 120 


VI  1  1 


. 


References  Cited  . . . 125 

7  .  APPENDIX  I  . ..132 

8.  APPENDIX  II  . 142 

9.  APPENDIX  III  . 166 


IX 


1: 


List  of  Tables 


Table 


Page 


1  Changes  in  Concentration  of  Organic  Carbon(OC): 

Comparison  of  Native  and  Cultivated  Sites . 14 

2  Changes  in  Mass  of  Organic  Carbon(OC): 

Comparison  of  Native  and  Cultivated  Sites . 17 

3  Effect  of  Crop  Rotations  on  Changes  in  Concentration 

of  Organic  Carbon(OC)  in  Soils . 21 

4  Effects  of  Additions  of  Organic  Residue  and 

Fertilizer  on  Concentration  of  Organic 

Carbon(OC)  in  Soil . 24 

5  Effect  of  Tillage  Intensity  on  Concentration  of 

Organic  Carbon(OC)  in  Soil . 26 


6  Significant  Correlation  Coefficients  Between 
Moisture  Retention  and  Management  and  Soil 
Variables  in  Cultivated  A  Horizons . 68 


7  Significant  Correlation  Coefficients  Between 

Moisture  Retention  and  Soil  Variables  in  Native 


A  Horizons . 69 

8  Relative  Frequencies  of  Management  Practices  on 

Farms  Sampled  in  Alberta  (  1982-1967) . 83 

9  Average  Annual  N  Inputs  (kg^ha-1)  to  Cereal 

Crops  Under  Various  Rotations  Between 

1982  and  1907 . 85 

10  Average  Annual  P  Inputs  (kg«ha_1)  to  Cereal 

Crops  Under  Various  Rotations  Between 

1982  and  1907 . 85 

11  Average  Annual  Yields  (Y)  (Mg^ha'1)  of  Cereal 

Crops  Under  Various  Rotations  Between 

1982  and  1907 . 87 


x 


12 


Average  Change  in  Concentration  and  Mass  of 
Organic  Carbon(OC)  in  A  Horizons  Under 
Summer  fallow ( SF )  and  Under  Continuous 
Cropping(CC) . 88 

13  Summary  of  Results  from  Multiple  Linear 

Regression  Models  Relating  Organic 
Carbon(OC)  Content  and  Change  to  Soil  and 
Management  Variables . 104 

14  Correlation  Coefficients  Between  Original 

Concentration  of  Organic  Carbon  and  Selected 
Management  Variables . 107 


xi 


■ 


F i gure 


List  of  Figures 


Page 


1  Three  Compartment  Model  of  Soil  Organic  Matter . 8 

2  Location  of  Sample  Sites  in  Alberta . 34 

3  Location  of  Sampling  Zones  in  Alberta . 34 

4  Schematic  Diagram  of  the  Method  of  Soil  Sampling . 36 

5  Relationship  Between  Inorganic  Carbon(%)  and  pH  of 

Soil  Samples  from  A,  B,  and  C  Horizons  of  Soils 
Throughout  Alberta . 40 

6  Change  in  Organic  Carbon  Based  on  Calculation  at 

Different  Levels  of  Hierarchy . 47 

7  Average  Concentration  of  Organic  Carbon  in 

Cultivated  and  Native  A  and/or  LFH  Horizons  of 
Soils  in  Alberta . 50 

8  Average  Mass  of  Organic  Carbon  in  Cultivated  and 

Native  A  and/or  LFH  Horizons  and'  Profiles  of 

Soils  in  Alberta . 52 

9  Average  Loss  in  Concentration  and  Mass  of  Organic 

Carbon  from  A  and/or  LFH  Horizons  of  Soils  in 
Alberta  Following  Cultivation . ,54 

10  Average  Thickness  of  Cultivated  and  Native  A 

and/or  LFH  Horizons  of  Soils  in  Alberta . 55 

11  Average  Apparent  Bulk  Density  of  Cultivated  and 

Native  A  and/or  LFH  Horizons  of  Soils  in  Alberta.. 57 

12  Average  pH  in  Water  of  Cultivated  and  Native  A 

and/or  LFH  Horizons  of  Soils  in  Alberta . . . 58 

13  Distribution  of  Organic  Carbon,  Apparent  Bulk 

Density  and  pH  with  Depth  in  Cultivated  and 

x  i  i 


Native  Soils  in  the  Brown  Sampling  Zone . 59 

14  Distribution  of  Organic  Carbon,  Apparent  Bulk 

Density  and  pH  with  Depth  in  Cultivated  and 
Native  Soils  in  the  Dark  Brown  Sampling  Zone . 60 

15  Distribution  of  Organic  Carbon,  Apparent  Bulk 

Density  and  pH  with  Depth  in  Cultivated  and 
Native  Soils  in  the  Black  Sampling  Zone . 61 

16  Distribution  of  Organic  Carbon,  Apparent  Bulk 

Density  and  pH  with  Depth  in  Cultivated  and 
Native  Soils  in  the  Dark  Gray  Sampling  Zone . 62 

17  Distribution  of  Organic  Carbon,  Apparent  Bulk 

Density  and  pH  with  Depth  in  Cultivated  and 
Native  Soils  in  the  Gray  Luvisolic  Sampling 
Zone . 63 

18  Distribution  of  Organic  Carbon,  Apparent  Bulk 

Density  and  pH  with  Depth  in  Cultivated  and 

Native  Soils  in  the  Peace  River  Sampling 

Zone . 64 

19  Average  Moisture  Retention  at  -33  and  -1500  kPa 

in  Cultivated  and  Native  A  and/or  LFH  Horizons 
of  Soils  in  Alberta . 66 

20  Average  Available  Moisture  Capacity  in  Cultivated 

and  Native  A  and/or  LFH  Horizons  of  Soils 
in  Alberta . 70 

21  Relationship  Between  Change  in  Apparent  Bulk 

Density  and  Change  in  Concentration  of  Organic 
Carbon  in  A  Horizons  of  Soils  in  Alberta . 76 

22  Relationship  Between  Change  in  Concentration  of 

Organic  Carbon  and  Time  in  the  Brown  Sampling 

Zone . 90 

23  Relationship  Between  Change  in  Concentration  of 

Organic  Carbon  and  Time  in  the  Black  Sampling 

Zone . 92 


XI  1  1 


24  The  Influence  of  Summer  fallowing  on  the 

Relationship  Between  Change  in  Concentration 
of  Organic  Carbon  and  Time  in  the 

Brown  Sampling  Zone . 94 

25  The  Influence  of  Continuous  Cropping  on  the 

Relationship  Between  Change  in  Concentration 
of  Organic  Carbon  and  Time  in  the 

Black  Sampling  Zone . 96 

26  Relationship  Between  Log  %  of  Original  Organic 

Carbon  and  Time  in  the  Brown  Sampling  Zone . 97 

27  Relationship  Between  Log  %  of  Original  Organic 

Carbon  and  Time  in  the  Black  Sampling  Zone . 98 

28  The  Influence  of  Summer  fallowing  on  the  Log  %  of 

Original  Organic  Carbon  and  Time  in  the  Brown 
Sampling  Zone . 99 

29  The  Influence  of  Continuous  Cropping  on  the  Log  % 

of  Original  Organic  Carbon  and  Time  in  the 

Black  Sampling  Zone . 100 


xiv 


1 .  Introduction 


Soil  organic  matter  is  a  continuum  of  materials  varying  in 
molecular  size,  reactivity  and  accessibility  (McGill , 1 983 ) . 
It  comprises  fresh  litter  and  plant  material  (tops  and 
roots)  in  various  stages  of  decay;  soil  biomass,  mostly 
fungi  and  bacteria  with  some  small  animals  and  insects  that 
decompose  residues;  and  microbial  by-products  of 
decomposition,  enzymes,  and  material  in  various  stages  of 
humification  (Campbell  and  Biederbeck,  1980). 

The  value  of  organic  matter  for  practical  agriculture 
flows  from  its  functions:  i)in  storing  and  supplying 
essential  plant  nutrients,  such  as  N,P,S  and  micronutrients 
through  mineralization  and  cation  exchange  reactions  (McGill 
and  Hoyt, 1977);  ii)in  maintaining  soil  structure  by 
influencing  aggregation;  and  thereby  water  infiltration, 
root  penetration,  seedling  emergence,  and  aeration 
( Schroeder , 1 978 ;  Scheffer  and  Schacht schabel , 1 979 ;  Lathwell 
and  Bouldin , 1 98 1 ) ;  iii)in  minimizing  erosion  via  surface 
residues  and  aggregation  (McGill , 1 982 ;  McGill  and 
Reinl,1984);  and  iv)possible  physiological  effects  (McGill 
and  Hoyt ,1977). 

These  functions  combine  to  influence  the  capability  of 
a  soil  for  crop  production;  and  upon  manipulation  by  man, 
its  productivity. 
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In  the  1920s  and  1930s,  researchers  started  to  document 
a  decline  in  organic  matter  on  the  Canadian  Prairies  with  a 
concomitant  decrease  in  yield  (Shutt,1925;  Caldwell  et 
al.,1939).  McGill  et  al.  (1981)  estimated  that  the  soils  of 
the  prairies  had  lost  35  -  50%  of  their  organic  matter,  this 
loss  having  three  fundamental  implications:  i)loss  in 
nutrient  supplying  power;  i i ) deter iorat ion  of  soil  physical 
properties;  and  i i i ) envi ronmental  problems  relating  to 
climate,  i.e.  an  increase  of  C02  levels  in  the  atmosphere. 

The  potential  economic  impact  of  the  loss  in  nutrient 
supplying  power  can  be  estimated.  For  example,  Beaton  (1980) 
reported  that  during  the  period  1974  to  1978,  grain  harvests 
contained  an  average  of  212,000  Mg  of  N  annually  in  Alberta. 
Fertilizer  addition  averaged  177,000  Mg  of  N  annually. 
Fertilizer  N  use  efficiency  seldom  exceeds  50%. 
Consequently,  N  removal  from  soils  averaged  at  least  124,000 
Mg  (212,000-177,000/2).  At  a  price  of  60  cents  per  kg  this 
equals  14.4  million  dollars  annually. 

Additional  concerns  arising  from  the  above  include  soil 
and  energy  conservation,  and  stability  and  profitability  of 
future  agricultural  production. 

The  best  estimates  now  available  suggest  serious 
organic  matter  depletion  problems  and  a  trend  toward 
continued  depletion.  They  are,  however,  based  upon  limited 
data  and  require  further  testing.  To  this  end,  a  Farming  for 
the  Future  project  supported  by  the  Alberta  Heritage  Trust 
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Fund  and  titled  "Conservation  of  Land  Productivity  in 
Alberta"  was  initiated. 

As  part  of  the  project,  this  thesis  has  the  following 
objectives : 

1 .  to  measure  the  extent  of  changes  in  organic  carbon 
in  cultivated  soils  of  all  climatic  zones  and  several 
cropping  systems  throughout  Alberta; 

2.  to  examine  interrelationships  among  observed 
changes  in  soil  organic  carbon  content ,  bulk  density , 
moisture  retention ,  and  regional  or  management 

f 

variables;  and 

3.  to  provide  a  broadly  based  data  source  for  testing 
existing  models  which  integrate  our  present  knowledge 
about  organic  matter  dynamics  and  which  will  be  useful 
for  future  management  decisions. 


2.  Literature  Review 


2 . 1  Introduction 


Sc ient i f ic 

interest 

in  soil 

organic 

matter 

has 

continued  since 

the  earliest 

sc ient  i  f  ic 

studies 

of  soils 

.  I  n 

1836,  the  Russian  scientist  German  found  that  cultivated 
soil  contained  less  humus  than  virgin  soil;  Dokuchaev  (1883) 
and  Kostychev  (1886)  reported  that  continuous  cultivation  of 
cereals  on  Russian  Chernozem  over  a  long  period  resulted  in 
decreased  soil  fertility  and  increased  susceptibility  to 
drought  (cit.  Kononova  et  al.,1966). 

These  findings  suggest  an  important  role  of  organic 
matter  in  crop  production  and  have  stimulated  many 
researchers  to  investigate  its  nature,  and  its  role  in  soil 
formation  and  fertility.  Past  efforts  have  concentrated  on 
chemical  and  physico-chemical  characterization  of  soil 
organic  matter;  with  interest  in  dynamics  developing  mainly 
in  the  last  two  decades.  The  dynamics  of  organic  matter 
require  kinetic  studies  and  descriptions  of  soil  systems. 
The  kinetic  approach  has  been 
questions  about  the  effect  of  cult 
and  its  change  with  time. 


useful  for  investigating 
ivation  on  organic  matter 
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2.2  Organic  Matter  Dynamics 

Organic  matter  normally  does  not  accumulate 
indefinitely  in  soil,  although  exceptions  can  occur,  as,  for 
example  in  permanently  waterlogged  soil.  Usually,  after  a 
period,  annual  losses  of  organic  matter  through 
decomposition  and  erosion  balance  annual  additions  and 
thereafter  the  organic  matter  content  of  the  soil  remains  at 
steady-state  ( Jenk inson , 1 966 ) .  This  process  in  which  losses 
and  gains  proceed  simultaneously  is  described  as  turnover, 
and  was  defined  as  flux  of  organic  material  through  the 
organic  matter  in  a  given  sample  of  soil  ( Jenk inson , 1 966 ) . 
Environmental  factors  which  influence  additions  and  removals 
and  thereby  the  balance  are:  temperature,  moisture, 
aeration,  and  pH  ( All i son , 1 973 ) .  Specific  combinations  of 
these  factors  determine  the  chemical  or  physical  nature  of 
organic  matter  formed. 

Changes  in  the  soil  environment,  such  as  breaking 

virgin  soil  or  altering  a  system  of  agriculture  cause 

changes  both  in  the  rate  of  addition  and  the  rate  of 

decomposition,  and  ultimately  a  new  steady-state  is  reached 
( Jenkinson , 1 966;  Lathwell  and  Bouldin , 1 98 1 ) .  The  direction 
of  change  depends  upon  the  antecedent  organic  matter  level 
as  well  as  the  nature  of  the  change  in  vegetation  system 
(McGill , 1 982 ;  Robert son , 1 98 3 ) . 


Man  influences  carbon  input  to  land  under  agricultural 
cultivation  by  controlling:  plant  species,  quantity  of  crop 
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grown,  proportion  of  time  during  which  a  crop  is  grown,  and 
crop  residue  management.  Output  is  influenced  by  controlling 
tillage,  soil  climate  and  erosion  (McGill  and  Reinl,1984). 

It  follows  that  productivity  of  cropping  systems 
directly  influences  the  level  of  organic  matter  at 
steady-state  (McGill  et  al.,1981).  Considering  the  various 
functions  of  organic  matter  in  soil,  it  also  follows  that 
management  practices  which  encourage  maintenance  of  soil 
organic  matter  should  aid  in:  improving  soil  structural 
stability;  improving  infiltration;  decreasing  evaporation; 
improving  drainage  in  fine  textured  soils;  and  fostering 
more  extensive  and  deeper  root  systems  which  in  turn 
increase  moisture  and  nutrient  availability  to  crops 
(Campbell  and  Biederbec k , 1 980 ) . 

Organic  matter  exists  in  forms  with  varying  degrees  of 
stability  and  accordingly  has  been  conventionally  divided 
into  three  components:  i)an  active  pool  with  a  turnover  time 
of  1-10  years  consisting  mainly  of  litter  and  biomass;  it  is 
the  pool  most  readily  affected  by  management  (McGill , 1 983 ) ; 
ii)a  pool  with  a  medium  turnover  time  of  about  60  years, 
termed  humads  by  McGill  and  Hoyt  (1977)  and  McGill  et 
al.(1981),  and  consisting  of  decomposition  products 
stabilized  by  adsorption  and  humification;  and  iii)a  pool 
with  a  turnover  time  of  a  thousand  years  or  more,  commonly 
referred  to  as  resistant  organic  matter. 

Several  mathematical  and  simulation  models  have  been 
used  to  describe  soil  organic  matter  dynamics.  They  have 
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been  reviewed  by  Van  Veen  et  al. ( 1981 )  and  in  Frissel  and 
Van  Veen(1981).  The  three  component  direct  input  model 
(McGill  and  Hoyt, 1977)  describes  organic  matter  dynamics  in 
terms  of  a  group  of  independent  first  order  reactions 
(Fig. 1 )  proceeding  concurrently  on  separate  substrates.  It 
has  been  used  to  provide  many  of  the  estimates  which  serve 
as  background  to  the  present  study.  Three  components  of 
different  decomposabi 1 i ty :  microbes,  humads  and  stable 
organic  matter  are  thereby  implicated.  Instead  of  flowing 
from  one  component  to  the  other,  carbon  flows  through  each 
component  independently. 

A  theoretical  objection  to  this  approach  is  the 
inherent  assumption  that  all  substrates  are  present  in  the 
residue  added.  Current  thinking,  however,  suggests  that 
considerable  transformation  of  one  substrate  into  another 
occurs  and  that  some  of  the  most  resistant  substrates  may 
not  be  stabilized  or  produced  until  some  time  after  the 
addition  of  residues  (McGill  and  Hoyt, 1977).  Restrictions  in 
use  and  interpretations  of  model  output  for  each  individual 
first  order  reaction  are:  i)addition  and  decomposition  rate 
must  remain  constant;  ii)turnover  time  of  the  total  soil 
organic  component  will  be  constant  only  in  a  steady-state 
situation;  and  iii)the  model  as  normally  used  suggests 
continuous  decline  of  organic  matter  with 
additions  are  very  high. 


time  unless 
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Fig.  1 


Three  Compartment  Model  of  Soil  Organic  Matter. 


9 


The  change  in  carbon  in  each  component  is  represented  by  the 
following  equations: 

dA/dt  =  -kaA  +  Ra  (  1  ) 

dB/dt  =  -kbB  +  Rb  (2) 

dC/dt  =  -kcC  +  Rc ;  Ra+Rb+Rc=1;  (McGill  and  Hoyt, 1977  ).  (3) 

Upon  integrating  ( Sr idharan , 1 970 ;  Tierney , 1 979 ) ,  the 
total  mass  of  carbon  at  time  t  (Mt)  is  determined  by  both 
decomposition  of  original  mater ial ( Ao , Bo , Co )  at  rates 
ka,kb,kc  and  by  additions  to  each  of  those  components  at 
rates  Ra , Rb  and  Rc : 

Mt  =  Aoe_kat  +  Ra/ka  (  1  -e ~kat  ) 

+  Boe~kbt  +  Rb/kb ( 1 -e ~kbt  ) 

+  Coe~kct  +  Rc/kc ( 1 -e ' kct  );  (McGill , 1 983  )  .  (4) 

A  different  quantity  of  organic  carbon  at  steady-state  is 
reached  for  each  component.  Steady-state  (Mss)  is  reached 
only  at  infinite  time,  i.e.  e~kat  approaches  zero,  but  the 
time  at  which  some  defined  proportion  (e.g.  0.9)  of 

steady-state  is  reached  can  be  calculated. 

Mt  =  0.9Mss  =  0.9(Ra/ka  +  Rb/kb  +  Rc/kc)  (5) 

if  Mt  is  increasing. 

For  component  A: 

0.9(Ra/ka)  =  Aoe~kat  +  Ra/ka  (  1  -e  "kat  ); 

0 . 9  ( Ra/ka  )-( Ra/ka  )(  1 -e "kat  )  =  Aoe_kat  ; 

Ra/ka  (  0 . 9- 1  (  1 -e " kat  )  =  Aoe'kat  ; 

Ra/ka  (  0 . 9- 1 +e  ~kat  )  =  Aoe"kat  ; 
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Ra/ka  (-0 .  1+e‘kat  )  =  Aoe"kat  ; 

-  O.IRa/ka  +  ( Ra/ka )  e  ~kat  =  Aoe'kat  ; 

-  O.IRa/ka  =  Aoe'kat  -  (Ra/ka)e'kat  ; 

-  O.IRa/ka  =  e_kat  (Ao  -  Ra/ka); 

-  0 . 1Ra/( ka (Ao-Ra/ka ) ) =  e‘kat 

-  0.1Ra/(Aoka  -  Ra )  =  e“kat  ; 

-  0 . 1 /( ( Aoka/Ra )  -  1)  =  e~kat  ; 

0  .  1  /  (  1  -(Aoka/Ra))  =  e'kat  ; 
ln(0.1/(1  -  (Aoka/Ra)))  =  -kat 

t  =  -1/ka  In ( 0 . 1/( 1 -Aoka/Ra )) ;  (6) 

The  same  holds  for  components  B  and  C.  The  time  needed 
to  reach  steady-state  is  a  function  of  the  decay  rate  (ka) 
and  of  the  difference  between  Ao  and  steady-state  (Ra/ka). 
The  decay  rate  dominates,  however,  because  of  the  In  Aoka/Ra 
term . 

% 

Consequently  component  A  reaches  its  steady-state  in 
5-10  years  whereas  component  C  requires  several  centuries. 
The  most  active  component  (A)  does  not  disappear  upon 
reducing  inputs  but  reaches  a  lower  steady-state  and  reaches 
it  earlier  than  do  components  B  or  C. 

Tracer  studies  have  shown  that  organic  materials  of 
very  different  chemical  constitutions  are  included  amongst 
those  most  readily  mineralized;  that  is  with  the  shortest 
turnover  time.  Such  isotopic  work  has  established,  however, 
that  physically  or  biologically  separable  fractions  of  soil 


organic  matter  are  involved  to  differing  extents  in  the 
turnover  of  organic  matter  (Greenland , 1 97 1 ) .  Consequently 
there  appear  to  be  important  physical  and  biological 
mechanisms  which  stabilize  organic  matter.  Clay  plays  an 
important  role  in  the  former  by  forming  organic  matter-clay 
complexes;  the  latter  involves  microorganisms  producing 
polysaccharides  in  pores  of  a  few  microns  diameter  enabling 
them  to  stabilize  microaggregates  effectively  from  within 
(Greenland, 1971 ;Foster, 1981 ) . 

Van  Veen  and  Paul  (1981)  concluded  from  a  computer 
simulation  of  organic  carbon  dynamics  in  grassland  soils 
that  the  steady-state  level  of  soil  organic  matter  was  more 
dependent  on  decay  rate  of  the  recalcitrant  fraction  than  on 
that  of  the  plant  residues  added  to  the  soil.  The 
steady-state  level(Mss)  is  further  dependent  on  the  addition 
rate  of  plant  residue  (R) ,  because  Mss=R/k;  a  change  in  R 
causes  a  proportional  change  in  Mss.  Van  Veen  and  Paul(1981) 
support  their  findings  with  a  statement  made  by  Jenny  in 
1930  (cit.  Van  Veen  and  Paul, 1981)  that  long-term  organic 
matter  levels  of  soil  are  primarily  dependent  on  the 
zonality  principle:  management  practices  of  cultivated  soils 
can  alter  soil  organic  matter  levels  as  long  as  the  practice 
is  maintained.  On  return  to  natural  conditions,  however, 
levels  will  reapproach  their  long-term  steady-state 
determined  by  the  soil  forming  factors  ( All i son , 1 973 ) .  Both 
decomposition  and  addition  are  influenced  by  climatic  zones 
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because  temperature  and  moisture  control  their  rates. 
Temperate  and  cool  humid  climates  favour  the  accumulation  of 
organic  matter,  whereas  hot  tropical  conditions  enhance  its 
decomposition  (Richards , 1 978 ) . 


In  situations  where  organic  matter  is  the 
of  -plant  nutrients,  emphasis  has  been  placed  on 
of  highly  active  forms  of  organic  matter.  With  i 
of  commercial  fertilizers,  as  is  especially 
intensive  agricultural  systems  of  Europe,  empha 
toward  the  well-humified,  resistant  fraction 
their  function  in  soil  aggregation,  exchange  ca 
as  a  source  of  minor  nutrients  ( All i son , 1 973 ) . 


main  source 
maintenance 
ncreased  use 
the  case  in 
sis  shifted 
s  because  of 
pacity,  and 


2.3  Effects  of  Cultivation  on  Organic  Matter 

Some  studies  on  changes  in  organic  matter  due  to 
cultivation  have  compared  cultivated  land  to  adjacent  native 
sites,  assuming  that  the  native  site  represents  original 


conditions.  Changes 

have 

been  reported  on 

either  a 

concentration  or  a 

mass 

basis;  seldom  both  at 

once.  Plot 

experiments  have  been 

used 

to  investigate  changes 

in  organic 

matter  as  a  func 

t  ion 

of  crop  rotations, 

f ert i 1 izer 

additions,  and  tillage 

inte 

nsity.  Tables  1  to  5 

summarize 

such  studies  with  an  emphasis  on  work  from  the  North 
American  Great  Plains. 
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Comparison  of  cultivated  and  native  sites  showed  that 
both  concentration  and  mass  of  organic  carbon  in  surface 
horizons  decline  after  land  is  cultivated  (Tables  1  and  2). 
Losses  in  concentration  range  between  14%  and  80%  depending 
on  soil  type,  management  and  cultivation  time  (Table  1). 
Losses  in  mass  range  between  12%  and  58%;  there  are  some 
gains  in  mass  reported  for  the  Luvisolic  and  Podzolic  Orders 
(Table  2).  Newton  et  al.(1945)  reported  that  the  average 
annual  loss  of  organic  carbon  from  older  cultivated  fields 
was  approximately  half  that  of  the  newer  cultivated  fields; 
the  rate  of  loss  decreased  with  increased  cultivation  time. 
The  same  observation  was  made  by  Shutt(1925)  and  in  other 
subsequent  studies. 

Brown  et  al.(1942)  found  significant  losses  also  in  the 
15-30cm  depth  except  for  Brown  Chernozemic  soils. 

Losses  were  attributed  to  decomposition  through 
microbial  activity  and  to  erosion.  The  extent  of  erosion, 
however,  is  only  speculative  in  most  cases.  Few  studies 
actually  measured  organic  matter  loss  due  to  erosion;  such 
work  was  reported  by  Slater  and  Carleton (  1 939 )  and  in 
succeeding  studies  conducted  at  the  Soil  Conservation 
Experiment  Station  at  Clarinda,  Iowa. 


assoc 

found 

soil 


Changes  in  bulk  density  have  been 
iation  with  changes  in  organic  matter, 
a  strong  inverse  relationship  between  bul 
organic  matter  (r=-0.908  between  bulk 


reported  in 
Russel ( 1960) 
k  density  and 
density  and 


Table  1  Changes  In  Concentration  of  Organic  Carbon(OC):  Comparison  of  Native  and  Cultivated  Sites. 
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nitrogen).  He  therefore  emphasized  the  importance  of 
determining  bulk  density  when  soil  changes  over  time  were 
considered,  particularly  when  changes  were  to  be  expressed 
on  a  mass  basis. 

Similarly,  Davidson  et  al.(1967)  concluded  from  a  study 
on  changes  in  organic  matter  and  bulk  density  with  depth 
under  2  different  cropping  systems  that  an  increase  in 
organic  matter  would  cause  a  decrease  in  bulk  density.  For 
an  organic  carbon  concentration  of  1.7%  the  authors  reported 
a  bulk  density  of  1.46  Mg  *  m~  3  ;  for  an  organic  carbon 
concentration  of  0.9%,  bulk  density  was  1.56  Mg-  irr  3  . 

Inverse  relationships  between  bulk  density  and  organic 
carbon  were  further  reported  by  Tanchandrphongs  and 
Dav  idson  (  1  970  )  ,  Black  (  1  973  )  ,'  Martel  and  Deschenes  (  1  976  )  , 
Jenkinson  and  Johnston ( 1 976 ) ,  Black  and  Siddoway ( 1 979 ) , 
Bauer  and  Black(1980),  Voroney  et  al.(1981),  Thiessen  et 
al.(1982),  and  Coote  and  Ramsey ( 1 983 ) . 

Management  effects  on  organic  matter  are  summarized  in 
Tables  3  to  5.  Rotations  containing  summerfallow  were  most 
destructive  of  organic  matter  (Shutt,1925;  Doughty  et 
al.,1954;  Hill, 1954;  Ridley  and  Hedlin,  1968;  Khan, 1969; 
Ferguson  and  Gorby,1971;  Pi ttman , 1 977 ;  Dormaar , 1 979 ;  Dormaar 
and  Pi ttman , 1 980 ) . 

Doughty  et  al. ( 1954)  demonstrated  the  destructiveness 
of  summer  fallowing  with  a  pot  experiment  using  samples  of  a 
cultivated  and  native  clay  loam  and  clay.  Organic  carbon 


21 


in 

o 

to 

c 


u 

o 


c 

0 

£ 

s_ 

cj 

o 

c 

ft 

o> 

£. 

o 


o 

c 

o 


io 

£_ 

c 

Q) 

U 

c 

o 

o 

c 


l/l 

0) 

O) 

c 

ro 

£ 

c_> 

c 

o 

in 

c 

o 


(0 

+-> 

o 

a 

a 

o 

£_ 

O 


o 

in 

•H 

a 

Q) 


4- 

UJ 


co 

a> 

n 

<a 


10 

in 

4-> 

0 

O 

a 

0 

Q. 

0 

c 

c 

0 

•  * 

o 

in 

o 

> 

0 

0 

**- 

4- 

< — 

0 

0 

4- 

4-’ 

4- 

4-» 

o 

*— 

•  * 

c 

•  - 

> 

T3 

o 

o 

0 

o 

0 

•  - 

0 

*0 

0 

c 

0 

r— 

ai 

*— 

4-* 

0 

0 

0 

c 

0 

4- 

•r- 

0 

0) 

> 

0) 

o 

d) 

o 

-H 

0)  V 

0 

4-< 

0 

+-* 

o 

o 

> 

*— 

0 

2: 

O) 

O) 

L. 

c 

0 

O 

•*- 

0 

0 

A 

i — 

-r~ 

L. 

•— 

o 

(0 

10 

0 

£_ 

•— 

in 

O) 

4-» 

4-» 

0 

0 

0 

c. 

CN 

c 

CM 

E 

0 

a 

0 

•*“ 

3 

0 

• — 

— 

0 

L. 

A 

d> 

d) 

E 

> 

u 

O 

t- 

0 

**- 

10 

c 

L. 

> 

T5 

> 

T3 

O 

< 

o 

u 

£_ 

l_ 

o 

o 

o 

**- 

10 

C 

10 

C 

o 

* 

-t— 

0 

'r“ 

CN 

in 

0 

in 

a 

' 

o 

— 

0 

0 

0 

c 

U 

l_ 

. — . 

L. 

o 

-—4. 

c 

-H 

in 

0 

C 

co 

0 

4-* 

O) 

0 

•H 

CN 

-H 

■D 

0 

co 

4-» 

n 

0 

CO 

L. 

3 

(J) 

< — 

C 

0 

O) 

0 

c 

• — 

O) 

0 

-C 

0 

0 

L. 

— 

< — 

<0 

t. 

4- 

00 

to 

C D 

— - 

in 

0 

0 

U 

cr 

in 

Cd 

T- 

r- 

>- 

co 

£. 

1 

1 

1 

1 

ai 

1 

c 

* 

• 

> 

1 

> 

0 

0 

a 

3 

o 

1 

£ 

-r- 

4- 

0 

•  - 

O 

■ — 

1 

-H 

r— 

i. 

c 

1 

"O 

a 

1 

o 

0 

0 

o 

•t- 

-H 

(0 

•  • 

i 

1 

•  - 

£ 

-H 

4- 

0 

c 

0 

a) 

c 

(0 

1 

C 

(0 

c 

0 

r— 

0 

L. 

c_ 

0 

E 

i_ 

1 

L- 

o 

4-> 

0 

£- 

0 

o 

o 

A 

o 

(0 

1 

o 

sz 

+ 

*•— 

a 

u 

c. 

u 

o 

1 

o 

3 

L. 

-M 

S- 

0 

0 

l 

d) 

1 

1 

1 

d) 

0 

> 

L. 

L- 

L. 

4-* 

C 

> 

-H 

c 

1 

•H 

c 

> 

4-* 

1 

> 

> 

> 

c 

L. 

o 

c 

£_ 

1 

c 

c_ 

O 

o 

o 

o 

O 

. — 

o 

o 

1 

o 

o 

- — 

q: 

CD 

CO 

o 

U 

u 

o 

o 

1 

1 

1 

o 

u 

O 

0 

CM 

in 

1 

O) 

• — ' 

in 

*•— 

CN 

to 

1 

o 

c 

c 

1 

■*- 

0 

-T— 

u 

+ 

+ 

co 

co 

1 

1 

1 

r 

c 

1 

+ 

1 

CJ 

o 

1 

o 

* 

* 

1 

1 

1 

1 

^ — - 

T3 

r- 

CO 

^r 

1 

in 

c 

CD 

o 

co 

1 

CD 

> — ' 

LU 

1 

o 

^r 

T- 

1 

o 

1 

1 

4- 

1 

o 

4-* 

*- 

* 

1 

1 

1 

o 

L. 

CN 

in 

in 

1 

CD 

c 

0 

CD 

in 

0) 

1 

O 

o 

4-* 

1 

u 

to 

co 

1 

1 

CM 

1 

1 

\ 

E 

1 

c 

u 

O 

in 

1 

in 

o 

in 

CO 

1 

N 

x: 

1 

r- 

1 

r- 

— 

4-* 

i 

in 

■*— 

1 

•»— 

a 

o 

T- 

i 

1 

i 

o 

0 

O 

1 

O 

X 

Q 

1 

1 

1 

r— 

1 

0 

1 

> 

> 

JZ 

1 

A 

t — 

T3 

a 

E 

D 

4-* 

E 

0 

> 

E 

1 

•— 

+J 

E 

— 

C 

■ — 

0 

r— 

0 

L. 

(0 

0 

1 

0 

< — 

<0 

o 

0 

0 

O 

u 

O 

0 

f— 

O 

1 

A 

— 

O 

00 

0 

o 

r“ 

0 

0 

2 

o 

1 

i n 

in 

* 

c 

(0 

1 

-r- 

o 

0 

0 

TJ 

1 

£_ 

> 

— 

-H 

0 

c 

1 

3 

c 

■H 

£_ 

1 

O) 

to 

— 

1 

O 

to 

<0 

0 

r 

~0 

0 

0 

c. 

1 

0 

£ 

o 

£) 

+■> 

-r- 

a 

10 

1 

0 

+J 

o 

•— 

0 

L. 

A 

X 

O 

»— 

1 

•*- 

d) 

< 

A 

o 

UJ 

»— « 

o 

1 

ca 

(/> 

in 

o 


c 

o 

V) 

c 

r 

o 

73 


o 


c 

L_ 

o 

u 


c 

o 

o 


CD 

I 


CD 

00 

CM 


03 
CO 

co  co  co 


CO 

r 

10  4-*  E 

c.  —  ro 
10—0 
£  in- 

(0 

T> 

C 

ro  c_ 

3  10 
O  — 

->  V 


.37  3.46  +  6  corn-oats-  et  a  1 .  significant 

clover  meadow  (1942)  only  under 

.39  3.36  -  1  cont .corn  cont.corn 


Location  Soil  Horizon/  Cone,  of  OC(°/o)  Change  Rotation  YRS  Reference  Comments 

Depth(cm)  Start  End  in 


22 


c 


—  in 

03 

£ 

•  2 

O) 

— - 

3 

03  2 

O 

O 

X 

2  o  x-  >  £ 

*->  O  +■'  O  fl>  (D  ■—  +-> 
to  x  c  —  in 

•m  C>ncOBBH)in 
O  -H  —  X  0)0  £  £  •-  —  O 

t_  O  <0  (0  C  +j4jJ3'- 

C  x  to  C 

X  C+'OCMt'UIO- 

>tJ  -  c  cr^r  w-t-'  tn  — 

■•-too—  OCDOOOO 
con  E  O't-t)’-  —  cl  o.  in 


L. 

c 

o 

to 

to 

CO 

4- 

0) 

JC 

X 

c 

£ 

4J 

d) 

o 

in 

3 

£ 

■r~ 

V 

X 

+J 

4-' 

•— 

■4-* 

0) 

o 

(0 

Q) 

C 

3 

4-» 

**- 

o 

o 

c 

O 

> 

o 

' — 

L. 

i — 

t. 

d) 

0) 

/ — « 

> 

*— 

o 

(0 

t> 

*— 

in 

CD 

c 

(0 

a> 

(0 

r 

C 

u 

(0 

(/> 

> 

0) 

CM 

0) 

I  I  I 


■*— 

T— 

•»— 

CO 

co 

CO 

03 

0) 

0) 

t> 

0) 

D) 


O 

c 

D 

O 

to 

0) 

£ 

4-* 

r— 

•  * 

o 

03 

03 

M-> 

0) 

E 

C 

c 

> 

■4-< 

ai 

-r- 

o 

-H 

O)  +■• 

4-» 

c 

(0 

c 

CO 

UJ 

0 

c 

£ 

* 

o 

0 

Q 

> 

c 

/'“'s 

0 

■*“ 

CO 

o 

0) 

r— 

TO 

■ — 

IS) 

L_ 

in 

“D 

c 

T3 

0) 

L_ 

0) 

•r» 

03 

ai 

▼— 

03 

•T- 

a: 

i 

w 

CD 

r» 

co 


co  o 
co  rr 


3 

O 

V 

<0 

tn  0)  tn  c 

■h  E  tn  x 

to  to  to  o 

0  X  4-  X  o 

I  Q)  1 —  O)  • 

C  >  <0  0) 

X  O  M-  3  C 

O  —  —  —  O 

O  O  to  C  O 


o 

c 

o 

CJ 


I 


o  --  cn 

■»“  04 


+  I  I 


CO 

in 

co 

CM 

co 

co 

CO 

co 

co 

CNJ 

tn 


i  +■» 


4— 

03 

a) 

0 

l/) 

O 

1 

L. 

C_ 

+-* 

1 

l 

2 

4- 

2 

D 

D 

03 

in 

tn 

4-  1 

»— 

I 

4— ’ 

4-' 

0 

4-* 

-H 

in  _c 

03 

1 

O 

in 

(/) 

r 

03 

03 

i  i 

1 

4- 

I 

03 

03 

2 

4- 

O 

O 

L.  4- 

O 

■ — 

4- 

2 

a 

a 

03 

l 

+ 

•  m 

l 

03 

in 

4-> 

4- 

1 

in 

u  i 

2 

1 

i 

4- 

4-' 

4— ' 

c 

0 

l 

03 

»  2 

l 

4- 

4- 

c  in 

c 

c 

0 

1 

1 

4- 

0) 

4-  1 

4- 

r— 

< — 

o  i 

o 

o 

o 

2 

in 

Cl 

in  ^ 

tn 

03 

03 

O  5 

o 

o 

in 

O 

co 

CD  CN 

o 

CN 

03 

co 

co 

in 

cn  in 

in 

m 

\ 

1 

i 

i 

1 

i 

1  1 

T— 

+ 

+ 

in 

LO 

^r 

CM 

in 

CM 

CN  CN 

o 

CN 

co 

CO 

CD 

-r_ 

‘T— 

rr  cn 

CN 

CM 

CN 

CN 

CM 

CN 

CO 

CO 

CD 

0) 

CO 

CO 

CO 

CO 

CM 

CO 

CO 

CO 

CO 

oi  co  tn  o  co  o 

•«-  <S>  O  04  co 


*  * 
00  CO 


*  * 
t"  T 

r-  co 


CO  OI  CO  CM  CO  CO 


tn  tn 


O  O 


in 

i 

o 


in 

t 

O 


in 


tn 

OI 


i 

O 


i 

o 


to 

r 

tn  *->  b 

x  —  to 

to  —  o 

2  in  — 


3 

C 

4-> 

tn  E 
to  to 
£  o 
o  — 


> 

0) 

o  > 
to  10 

co  o 


a 

o 

—  -o 
o  c 
O  to 
c_>  tn 


to 

t> 

c 


to  X 
3  to 
o  — 
u 


to 

a; 

£_  I  O) 

mro 

.0  -H  — 
—  0)  X 
<  _i  XJ 


(0  O) 

xi  ai 
o  Q. 

■4-*  T- 

—  c 
c  c 
(0  — 
2  3 


C  1 
£-  to 

0)  X 


in  tn  to 
a>  3  — 

3  <  - 


. 


23 


TJ 

0 

3 

C 


C 

O 

u 


co 

0) 

n 

0 


c 

0) 

E 

E 

O 

c_> 


0) 

u 

c 

d) 

c_ 

0 

M- 

0 

cc 


to 

a. 

>■ 


to 


\  E 
c  u 

O  w 

n  sz 
■—  -*-> 
t_  a 
o  0 

I  Q 


O 

to 


C 

o 


(0 

o 

o 


C  0) 

0  u> 
r  m 
m  — 


a> 

CO 


c 

0 


a 

< 


o 

to 

>- 
0  > 
£-  3 
a  _i 


to 

•H  c 
C-  o 

0  -H 

23  0 

—  C_ 
<  CD 


>4-  to 

o  c 
o 

0  — 
U)  -M 
0  0 
C-  -t-1 
0  O 
>  c. 
0 
* 


C  —4 
0  r~- 
E  t" 
•H  CD 


O 

CM 


(D 

E 

1 

E 

+ 

1 

U) 

D 

1 

c 

(0 

O) 

0) 

1 

0 

<D 

c 

<D 

QJ 

c 

• 

0 

E 

£_ 

o 

r~ 

•— 

1 

r 

+ 

+ 

0) 

4- 

(D 

— 

i 

l 

E 

1 

4- 

3 

4- 

y 

W 

4- 

+J 

C 

C 

+ 

1 

(/) 

W 

i 

ID 

r— 

0 

4- 

**- 

4- 

1 

1 

1 

y 

5- 

4- 

-H 

<0 

■H 

(0 

(fl 

<D 

(D 

1 

5 

c 

i 

> 

r— 

c 

4- 

o 

L. 

1 

s_ 

1 

1 

1 

o 

1 

4- 

o 

■ — 

o: 

0) 

5 

0) 

1 

2 

a 

w 

CD 

+ 

u 

ro 

1 

0 

s? 

■* 

* 

1 

O) 

w 

1 

O 

CM 

CM 

CN 

CD 

CD 

c 

c  • 

T— 

CN 

1 

^r— 

— ■ 

T— 

*— 

CN 

T— 

0 

—  u 

1 

1 

1 

1 

i 

1 

i 

1 

1 

o 

c 

1 

o 

o 

1 

o 

1 

1 

1 

1 

-X- 

* 

"O 

^r 

▼“ 

r- 

CD 

1 

CD 

05 

CD 

O 

co 

CD 

c 

1 

CD 

■»T 

in 

r- 

O 

ID 

1 

*»— 

CN 

o 

1 

*»“ 

— 

•*- 

■*- 

in  O 

4“  co 
i  i 

o  o 


c 

*  3 

l-  o 
0  c 

O  CD 


CO  05  O  — 
CD  CD  t"  r~ 


in 


0 

4-> 

£.  I 
0  SZ 
£>  +-> 
—  0 
<  _l 


0 

O) 

T3 

£_ 

12 


C_ 

o 

-H 

0 

Ifl  uu 
0  —  0 
•a 

z  C  <4- 
*  —  o 


c 

o 

10 

3  T2 
ai  c 

t-  0 

0 


CN  00 


>  — 
n  c~ 
t_  cn 
o  — 

o  w 


o  o 


*  * 
cm  in 
co  or 

CM  CM 

O  6 


in 

i 

O 


n 

0 


0 

> 

0 

+-" 

o 

c 

c 

o 


0 

E 

£_ 

o 


05 


M- 

c 

in 

in 

i 

c 

>  0 

3 

3 

3 

o 

• 

— 

•— 

O 

0 

O 

c 

o 

£_ 

— 

3 

• — 

0 

L. 

— 

Q) 

w 

c 

C 

■ — 

C_ 

0 

E 

■o 

c 

c 

o 

(D 

— 

0 

0 

0 

sz 

0 

L_ 

E 

0 

4- 

4-* 

£- 

<4- 

u 

o 

N 

o 

3 

■O 

«— 

c 

3 

L. 

• — 

■ — 

(0 

o 

C 

0 

0 

0 

O 

JC 

4-> 

4- 

o 

0 

E 

D 

a 

C 

u 

o 

• — 

c 

E 

in 

E 

0 

o 

o 

3 

C- 

ID 

L. 

>« 

4-* 

0 

3 

0 

*- 

T5 

£1 

o 

0 

<T3 

>  V) 

SZ 

C 

c 

•  • 

a 

V 

r 

- 

O 

t- 

3 

0 

c 

0 

0 

•  • 

in 

4- 

c 

C 

*— 

0 

£- 

. — 

23 

a 

»— 

o 

-• 

0 

•• 

4- 

0 

s 

CD 

CQ 

O 

> 

(0 

a 

u 

2: 

E 

0 

L. 

in 

3 

0 

1 


Table  4  Effects  of  Additions  of  Organic  Residue  and  Inorganic  Fertilizer  on  Concentration  of  Organic  Carbon(OC) 
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content  was  measured  at  the  beginning  and  end  of  5  years. 
Losses  varied  from  18.6%  for  the  cultivated  Sceptre  clay 
(1.94%  start;  1.58%  end),  19.9%  for  the  cultivated 
Woodmountain  clay  loam  (2.73%  start;  2.26%  end)  to  21%  for 
the  virgin  clay  (2.73%  start;  2.15%  end)  and  24%  for  the 
virgin  clay  loam  (3.8%  start;  2.88%  end).  The  higher  losses 
in  the  virgin  soils  were  attributed  to  a  higher  proportion 
of  more  easily  decomposable  material. 

Ridley  and  Hedlin(1968)  implied  that  rapid  decline  of 
organic  matter  in  frequently  fallowed  land  was  due  to  higher 
frequency  of  tillage  operations.  Tillage  was  considered  to 
disrupt  aggregates  and  expose  new  areas  to  microorganisms 
thus  promoting  more  rapid  oxidation  of  organic  carbon  as 
suggested  by  Rovira  and  Graecen ( 1 957 ) . 

More  intensive  cropping  slowed  down  the  rate  of  organic 
matter  decline  which  was  attributed  to  greater  addition  of 
residue  in  the  form  of  roots  because  aboveground  residues 
had  been  removed  (Ridley  and  Hedl in , 1 968 ) .  Roots  are 
intimately  mixed  with  the  soil  at  all  times.  Their 
excretions  and  residues  provide  substrate  for  microorganisms 
(Allison , 1 973 ) .  Extracellular  microbial  products  participate 
in  stabilization  of  microaggregates  and  slow  decomposition 
by  physical  protection  (Greenland , 1 97 1 ) . 

Grasses  and  legumes  decreased  losses  in  organic  matter 
(Blair  and  McLean, 1917;  Shutt,1925;  Salter  and  Green, 1933; 


I 
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Johnston  et  al.,1942;  Van  Bavel  and  Schaller , 1 950 ; 
Hill, 1954;  Khan, 1969;  Sowden  and  Atkinson , 1 968 )  through  high 
below  ground  primary  production  (McGill  and  Hoyt, 1977). 

Jenkinson ( 1 977 )  monitored  the  rate  of  carbon  loss  over 
10  years  from  C14-labelled  plant  material  incubated  in  soils 
under  fallow  and  grass.  Decomposition  under  growing  grass 
was  considerably  slower  than  under  bare  fallow.  After  5 
years  there  was  47%  more  of  the  originally  added  material 
under  grass  than  under  fallow.  The  protection  given  by  the 
grass  ended  as  soon  as  it  was  removed  leading  to  a 
decomposition  rate  of  the  labelled  material  which  was  twice 
as  high  as  under  bare  fallow  and  after  10  years  both  soils 
contained  approximately  the  same  amount  of  added  material. 
The  slower  decomposition  rate  under  grass  could  not  be 
explained  by  recycling  of  C14  through  photosynthetic 
fixation,  translocation,  and  return  to  the  soil  because  C14 
content  of  the  grass  was  too  low  to  be  detected. 

The  rapid  accumulation  of  soil  organic  matter  usually 
observed  when  arable  land  is  put  under  grass  is  attributed 
to:  i) increased  annual  return  of  carbonaceous  material;  and 
ii)decreased  rate  of  decomposition  under  grass 
(Jenkinson , 1 977 ) .  The  latter  was  ascribed  to  limiting  N 
supply  (Huntjens  and  Albers , 1 978 ) . 


Row 

organic 

Johnston 


crops  appeared  to  be  more  destructive  to  soil 
matter  than  cereal  crops  (Salter  and  Green, 1933; 
et  al.,1942;  Coote  and  Ramsey , 1 983 ) .  Row  crops 
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require  frequent  tillage  for  reducing  weed 
soil  is  left  bare  during  a  longer  portion 
period  due  to  slow  initial  growth 
(Geisler , 1 980 ) . 


competition;  and 
of  the  vegetation 
of  row  crops 


Additions  of  straw  and  manure  may  increase  organic 
matter  (Table  4).  Manure  is  a  source  of  the  major  plant 
nutrients  N,P,K  (Allison , 1 973 )  and  like  other  residues  may 
ameliorate  soil  conditions  for  plant  growth  thus  enhancing 
primary  productivity  (McGill  and  Hoyt, 1977). 

Application  of  mineral  N  and  P  seems  to  be  conducive  to 
preventing  organic  matter  losses  by  increasing  production 
and  therefore  return  of  residues,  mainly  of  roots  (Jenkinson 
and  Johnston , 1 976 ) . 


Zero  till,  minimum  till  and  stubble  mulch 
the  rate  of  organic  matter  loss  as  opposed  to  c 
tillage  practices  (Table  5).  Aeration,  oxidation 
are  reduced  because  little  physical  disruption  of 
occurs.  Doran ' s ( 1 980 )  data  indicate  that  the  posi 
of  zero  tillage  might  be  mainly  present  in  surfac 
There  may  be  a  greater  loss  in  subsurface  hor 
zero  tillage  than  under  conventional  tillage, 
attributed  this  to  placement  of  crop  residues  at 
ploughing  under  conventional  tillage. 


ing  lowered 
onvent ional 
and  erosion 
aggregates 
tive  effect 
e  horizons, 
izons  under 
The  author 
depth  with 
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2.4  Conclusions 

Magnitude  and  direction  of  change  in  organic  matter  are 
related  to  the  balance  between  original  amount  of  organic 
matter  and  kind  of  management  introduced. 

The  rate  of  organic  matter  loss  tends  to  be  highest  in 
the  initial  years  of  cultivation;  thereafter  losses  continue 
at  lower  rates.  Losses  in  mass  tend  to  be  lower  than  losses 
in  concentration  which  has  been  reported  to  be  due  to  an 
inverse  relationship  between  organic  carbon  and  bulk 
density . 

Losses  occur  if  the  removal  rate  exceeds  the  addition 
rate.  Frequent  tillage  operations  and  summer  fallowing 
increase  removal  rates  directly  by  enhancing  decomposition 
rates;  indirectly  by  making  the  soil  more  susceptible  to 
erosion.  Decomposition  and  erosion  are  linked  by  a  positive 
feedback  mechanism. 

Rotational  studies  showed  that  organic  matter  may 
increase  if  legumes  and  grasses  are  grown;  and  that  it 
decreases  most  rapidly  if  summer  fallowing  is  practiced. 
Continuous  cropping  slows  organic  matter  losses  mainly  by 
increased  addition  of  root  residues.  Additions  of  organic 
residues  increase  the  addition  rate  and  may  increase  organic 
matter  if  the  removal  rate  remains  constant.  Zero  till  is 
considered  to  reduce  decay  rates  and  thus  cause  organic 
matter  levels  to  rise.  Data  of  Doran(1980)  indicate  that 
this  might  be  a  redistribution  effect  rather  than  actual 
accumulation  of  organic  matter. 
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It  follows  from  the  above  that  a  change  in  management 
may  cause  organic  matter  losses  to  cease  or  change  direction 
if  the  introduced  management  practice  either  lowers  removal 
rate,  e.g.  by  reduced  tillage  intensity,  or  increases 
addition  rate. 

Very  few  such  measurements  have  been  made  on  farm 
fields  so  far.  Mathematical  treatments  of  organic  matter 
dynamics  under  the  influence  of  cultivation  are  based  on  the 
assumption  that  management  remains  constant.  Therefore,  in 
their  present  form,  they  are  inadequate  to  describe  systems 
where  changes  in  management  alter  organic  matter  dynamics. 


3.  Materials  and  Methods 


3 . 1  Field  Work 


Six  major 
delineated 


Soil  samples  were  collected  from  cultivated  fields  and 
adjacent  native  sites  on  72  farms  located  in  various 
climatic  and  soil  zones  of  Alberta  (Fig. 2) 
sampling  zones,  from  south  to  north,  were 
(Fig . 3 )  : 

1 .  Brown  ( BR) 

2.  Dark  Brown  (DB) 

3.  Black  (BL) 

4.  Dark  Gray  (DG) 

5.  Gray  Luvisol  (GL)  soil  zones  and  the 

6.  Peace  River  area  (PR)  which  contains  soils  of  the  Solod, 
Solodized  Solonetz,  Gray  Luvisol  and  Dark  Gray  Great  Groups. 


3.1.1  Site  Selection 

Farms  were  selected  with  the  help  of  Alberta 
Agriculture  through  District  Agriculturalists,  Regional 
Plant  Industry  personnel  and  Regional  Soils  Specialists. 
They  contacted  and  provided  names  of  farmers  who  were 
willing  to  cooperate  and  who  possessed  possible  sample 
sites . 

The  following  criteria  were  applied  in  finally 

selecting  sampling  sites: 

1.  location  of  a  cultivated  field  adjacent  to  an  undisturbed 
native  site  accessible  to  a  coring  truck; 
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Fig .  2 
Fig .  3 


Location  of 
Location  of 


Sample  Sites  in  Alberta. 
Sampling  Zones  in  Alberta. 
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Fig. 2  Location  of  Sample  Sites  in  Alberta. 
Fig. 3  Location  of  Sampling  Zones  in  Alberta 
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2.  presence  of  soils  of  the  same  origin  on  both  sites  and 
therefore  of  the  same  parent  material,  microclimate, 
topography,  drainage,  and  native  vegetation; 

3.  availability  of  field  history  and  good  management 
records ; 

4.  absence  of  effects  due  to  fencelines,  wind-drifts,  and 
former  threshing  piles. 


The  distribution  of  sites 
uneven  due  to  different  farmin 
sites  within  sampling  zones  wa 
and  adequate  sampling  condi 
density,  the  site  numbers  had 
and  financial  restrictions, 
the  geographic  distribution  wa 
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3.1.2  Soil  Sampling 

At  each  farm,  8 
hydraulic  soil  corer 
The  diameter  of  the  c 
depth  was  one  meter, 
possible.  Problems  ar 
sufficient  depth  wa 
consisted  of  a  lar 
compaction.  Apparent 
to  the  top  of  the  soi 


prof i le 
along  a  t 
oring  tube 
Sampling 
ose  on  roc 
s  difficul 
ge  propor 
compaction 
1  in  the  c 


pairs  were  sampled  with  a 
ransect  as  depicted  in  Fig. 4. 

was  4.5  cm,  and  its  maximum 
was  to  the  C  horizon  wherever 
ky  soil  as  penetration  to 
t  or  when  subsurface  material 
tion  of  clay  which  caused 
from  the  surface  of  the  soil 
oring  tube  was  measured  for 
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O  Soil  Core  -  Lab  Analysis 
•  Soil  Core  -  Density  Sample 


Cultivated 


Fig .  4 


Schematic  Diagram  of  the  Method  of  Soil  Sampling. 


For  five  of  the  profile  pairs, 
delineated,  horizon  thickness  measured, 


horizons 
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were 


texture  of  the 


parent  material  determined  and  each  horizon  sample  stored 
separately  for  analyses. 

At  15  to  30  cm  away  from  cores  #1,#3,  and  #5,  an 
additional  core,  was  taken  for  the  purpose  of  determining 
mass  and  apparent  bulk  density  for  each  horizon.  Horizons 
were  therefore  left  undisturbed;  compaction  and  horizon 
depth  were  measured  and  each  horizon  sample  stored 
separately . 


3.1.3  Field  History 

A  questionnaire  (Appendix  I)  was  provided  to  each 
farmer  and,  whenever  possible,  completed  on  site  with  the 
farmer.  The  information  derived  consists  of  the  total  time 
of  cultivation,  crop  rotations,  management  practices,  and 
yields,  aiming  at  the  assessment  of  input-output 
relationships  for  the  respective  soil-cropping  system. 

Questionnaires  which  could  not  be  answered  in  the  field 
were  subsequently  mailed  in.  Some  questionnaires  were  not 
returned  and  a  complete  statistical  analysis  could  not  be 
performed  in  those  cases. 


/ 
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3.2  Lab  Procedures 

3.2.1  Sample  Preparation 

Samples  were  air-dried  and  their  mass  determined  on  a 
top-loading  balance.  Samples  for  lab  analysis  were  ground  to 
2  mm.  A  subsample  was  ground  in  a  Siebtechnik  laboratory 
disc  mill  (Model  TS100A)  to  100  mesh  which  is  equivalent  to 
a  particle  size  of  0.149  mm  or  very  fine  sand  according  to 
the  Wentworth  size  class  (Quaternary  Research  Group, 1979). 

For  each  sample  set  the  dry  Munsell  color  of  the 
A-horizon  was  recorded  for  classification  purposes. 

3.2.2  Apparent  Bulk  Density 

Bulk  density  is  the  ratio  of  mass  to  macroscopic  volume 
of  soil  particles  plus  pore  spaces  in  a  sample  (Black  et 
al . , 1 965 , Vol . 1 ) .  The  mass  of  the  samples  taken  for  bulk 
density  measurement  was  therefore  divided  by  the  product  of 
horizon  depth  and  area  of  the  coring  tube  (15.9 
cm2 ). Compact  ion  measured  was  minor  except  for  a  few  C 
horizons  and  was  not  included  in  the  calculation  of  bulk 
density.  Because  of  this,  together  with  the  fact  that 
air-dry  instead  of  oven-dry  mass  was  used  in  the 
calculation,  bulk  density  will  be  referred  to  as  apparent 
bulk  density. 
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3.2.3  Soil  pH 

Soil  pH  was  measured  in  water  for  the  2  mm  samples  as 
described  in  McKeague  ( 1 978 ; p67-68 ) .  The  suspension  had  a 
soilrwater  ratio  of  1:2.5;  this  value  was  adopted  in  1930  by 
the  Soil  Reaction  Committee  of  the  International  Society  of 
Soil  Science  (Black  et  al.,1965). 

Measurements  were  conducted  potent iometr ically  with  a 
Fisher  Acumet  pH  meter  Model  630.  The  instrument  was 
standardized  in  buffer  solutions  at  pH  4  and  pH  7.  The  glass 
electrode  was  immersed  into  the  partly  settled  solution,  the 
calomel  electrode  into  the  supernatant  to  reduce  errors 
caused  by  the  liquid  junction  potential  in  clay  and  soil 
suspensions  (Black  et  al.,1965).  Results  were  recorded  to 
one  decimal. 

3.2.4  Inorganic  Carbon 

Inorganic  carbon  was  determined  for  2  mm  materials  by 
acid-neutralization  in  boiled  samples  (Black  et 
al . , 1 965 ; pi  387- 1  388 )  whenever  the  pH  of  a  sample  was  equal 
to  or  greater  than  7  which  was  the  case  for  800  samples. 
Plotting  of  the  carbonate  values  versus  the  pH  measurements 
(Fig. 5)  indicates  that  a  pH  of  7  is  acceptable  as  a  cut-off 
point  for  carbonate  determination  for  purposes  of  the 
present  study. 

The  recommended  molarity  of  the  NaOH  (0.25  M)  was 
inconvenient  for  samples  with  low  carbonate  content  because 
of  the  high  amount  of  base  needed  to  neutralize  the  acid 
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pH 


Fig. 5  Relationship  Between  Inorganic  Carbon(%)  and  pH  of 
Soil  Samples  from  A,  B,  and  C  Horizons  of  Soils 
throughout  Alberta  (n=800). 
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(HC1).  The  concentration  of  NaOH  was  therefore  raised 
to  0.5  M.  The  actual  concentration  of  base  was  determined 
daily  by  titrating  a  known  volume  of  base  with  standard  HC1, 
in  order  to  overcome  changes  in  NaOH  concentration  caused  by 
absorption  of  H20  and  C02 .  If  samples  showed  violent 

frothing  with  HC1  due  to  high  carbonate  content,  further 
measured  amounts  of  HC1  were  added  until  no  sign  of  a 
reaction  occurred. 

Visual  recognition  of  the  phenolphthale i n  endpoint  was 
difficult  due  to  the  often  yellow-brown  color  of  the 

filtrates  which  were  collected  after  the  soil-acid  mixture 
had  been  boiled  and  cooled.  The  method  was  therefore  changed 
by  using  a  potent iometric  titration  to  pH  8.2  controlled 
automatically  by  a  Radiometer  Titrator  Type  TTT  lib  operated 
according  to  the  manufacturer's  manual.  The  two  methods 
yielded  the  same  results  when  compared  on  the  same  samples. 

A  blank  was  run  with  each  set  of  samples.  The 

coefficient  of  variation  of  an  internal  standard  soil  sample 
was  2.26%. 

3.2.5  Total  Carbon 

Total  carbon  was  determined  by  oxidation  at  1300°C  with 
a  resistance  furnace  (Leco  Carbon  Determinator  CR12,  Model 
780-000)  in  an  02  stream  followed  by  infra-red  C02  detection 
as  described  in  the  manufacturer's  manual.  Mass  of  the 
sample  (<100  mesh  size)  varied  between  0 . 1 g  and  0.3g 
depending  on  the  expected  carbon  content.  Previously  dried 
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CaC03(12%)  stored  in  a  dessicator  was  used  for  calibration. 
The  coefficient  of  variation  of  an  internal  standard  soil 
sample  was  3.5%.  Results  were  recorded  to  two  decimals. 

3.2.6  Organic  Carbon 

Organic  carbon  was  calculated  by  subtracting 
%  inorganic  carbon  from  %  total  carbon  and  reported  as  such 
rather  than  converted  to  organic  matter  content  because 
considerable  variation  exists  in  the  Crorganic  matter  ratio 
among  different  soils,  and  also  among  horizons  of  the  same 
soil  (Black  et  al ., 1 965 ; Page  et  al.,1982). 

3.2.7  Moisture  Retention 

Moisture  retention  at  ~33kPa  and  -1500kPa  was 
determined  for  disturbed  2  mm  samples  of  A-horizons.  A 
Pressure  Plate  Extractor  CAT. #1600,  Soil  Moisture  Equipment 
Co., Santa  Barbara,  California  was  used  for  -33kPa  and  a 
Ceramic  Plate  Extractor  CAT. #1500,  Soil  Moisture  Equipment 
Co., Santa  Barbara,  California  for  -1500kPa,  respectively. 

The  sample  size  was  25g,  retainer  rings  had  a  diameter 
of  5cm.  A  control  soil  was  run  with  each  sample  set.  The 
coefficient  of  variation  of  the  control  at  -33kPa  was  11.7% 
and  18.4%  at  -1500kPa. 
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3.3  Data  Management 

Field,  lab  and  questionnaire  data  were  coded  and 
recorded  on  coding  sheets  prior  to  keypunching,  entry  into 
three  MTS  files,  and  further  preparation  for  statistical 
analysis.  The  contents  of  the  files  were  respectively: 

i)  farm  location  and  pedological  site  description; 

ii)  management  and  input-output  data  in  terms  of  yields  and 
fertilizer  additions;  and  iii)  chemical  and  physical 
analyses  data.  Data  were  stored  on  tape.  The  codebooks  are 
listed  in  Appendix  I.  Statistical  analyses  used  SPSS  and 
SPSSx  statistical  packages. 


4.  Changes  in  Soil  Properties 


4 . 1  Introduction 

The  objective  of  this  section  is  to  document  changes  in 
soil  properties  due  to  cultivation.  The  following  soil 
characteristics  are  compared  on  cultivated  and  native 
samples:  concentration  and  mass  of  organic  carbon  and  their 
respective  changes;  horizon  thickness;  pH;  apparent  bulk 
density;  moisture  retention  at  -33  and  -ISOOkPa;  and 
available  moisture  capacity.  The  change  in  organic  carbon 
when  calculated  at  the  profile,  farm  and  sampling  zone  level 
is  also  presented. 


4.2  Methods  of  Statistical  Analysis 

Data  presented  and  discussed  in  this  section  were 
obtained  from  lab  analyses  and  field  measurements.  Using 
SPSSx  subprogram  T-Test  the  null  hypotheses  of  no  difference 
between  the  following  means  have  been  tested: 

1.  Concentration  of  organic  carbon  in  cultivated  and  native 
A  horizons. 

2.  Mass  of  organic  carbon  in  cultivated  and  native  A 
hor i zons . 

3.  Mass  of  organic  carbon  in  cultivated  and  native 

profiles(to  bottom  of  B  horizon). 

4.  Change  in  organic  carbon  on  a  concentration  and  mass 

basis. 

5.  Depth  of  cultivated  and  native  A  horizons. 
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6.  Apparent  bulk  density  of  cultivated  and  native  A 
horizons . 

7.  The  pH  of  cultivated  and  native  A  horizons. 

8.  Distribution  of  organic  carbon,  apparent  bulk  density, 
and  pH  with  depth  in  cultivated  and  native  profiles. 

9.  Moisture  retention  at  -33  and  -1500kPa,  and  available 
moisture  capacity  in  cultivated  and  native  A  horizons. 


At  the  farm  level  t-tests  were  based  on  the  assumption 
that  samples  were  independent.  Pooled  or  separate  sample 
variances  were  used  depending  on  whether  samples  came  from 
populations  with  equal  or  unequal  variances  determined  by  an 
F~Test(Nie  et  al . , 1 975 ,p. 270 ) . 

A  paired  t-test  was  used  for  testing  null  hypotheses  on 
a  soil  zone  level.  Because  of  the  influence  of  soil  forming 
factors  it  was  assumed  that  means  of  cultivated  and  native 
sites  per  farm  would  be  more  closely  related  to  each  other 
than  would  be  means  of  cultivated  and  native  sites  of  all 
farms  within  a  zone.  Thus,  small  but  consistent  differences 
were  more  likely  to  become  statistically  significant. 

Highly  significant,  and  significant  are  used  throughout 
to  designate  statistical  significance  at  the  1  %  and  5  % 
level , respectively. 


Results  are  presented  for  the  sampling  zone  level  only. 
Summary  results  by  farm  are  included  in  Appendix  II.  In  each 
of  the  Brown(BR),  Dark  Brown(DB),  Black(BL),  Dark  Gray(DG), 
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Gray  Luvi sol ic (GL )  and  Peace  River(PR)  sampling  zones  14, 
16,  15,  2,  7,  and  12  farms  were  analyzed  statistically.  This 
amounts  to  70,  80,  75,  10,  35  and  60  profiles  at  each  of  the 
cultivated  and  native  sites  in  the  above  sampling  zones, 
respectively . 


4.3  Results 

Content  of  organic  carbon(OC)  can  be  expressed  as 
concentration  (g/IOOg  soil)  or  mass  (Mg.ha"1).  The  mass  of 
organic  carbon  per  hectare  over  a  given  depth,  e.g.  depth  of 
the  A  horizon,  is  obtained  by  multiplying  the  concentration 
by  bulk  density  and  horizon  depth,  accounting  for  changes  in 
units.  Concentration  indicates  the  amount  of  organic  carbon 
associated  with  mineral  soil  constituents  which  is  of 
interest  with  regard  to  soil  physical  properties  such  as 
aggregation;  it  has  a  qualitative  connotation.  Mass  is 
pertinent  in  terms  of  quantitative  aspects  for  assessing 
nutrient  release,  release  of  C02  to  the  atmosphere,  or 
organic  carbon  lost  by  erosion. 


4.3.1  Calculation  of  Changes  in  Organic  Carbon 

The  possible  range  of  results  obtained  when  calculating 
changes  in  organic  carbon  at  different  levels  of  hierarchy 
such  as  profile,  farm,  and  soil  zone  is  illustrated  in 
Fig. 6.  An  example  calculation  for  the  GL  zone  is  presented 
in  AppendixII.5  (p.158).  The  presence  of  a  denominator  in 
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Concentration  Q  Zone 
YA  Farm 

EQ  Individual  Core 
Mass  □  Zone 

E3  Farm 

ED  Individual  Core 


Fig. 6  Change  in  Organic  Carbon  Based  on  Calculation  at 
Different  Levels  of  Hierarchy. 
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the  equation  for  organic  carbon  change1  can  skew  the  results 
depending  on  the  level  of  hierarchy  and  variability  in  the 
denominator : 

1.  Changes  can  be  calculated  directly  from  the  difference 
between  the  means  of  organic  carbon  at  the  sampling  zone 
level.  These  means  are  derived  from  the  average  organic 
carbon  content  of  cultivated  and  native  sites  at  the  farm 
level.  This  method  does  not  provide  for  calculation  of 
standard  deviations  of  the  change. 

2.  Changes  can  be  calculated  at  the  farm  level  and  averaged 
for  the  soil  zone  level. 

3.  Changes  can  be  calculated  on  the  profile  level  and 
averaged  for  the  farm  and  soil  zone  level. 

Normally  changes  are  most  negative  with  the  first 
method  and  least  negative  with  the  third.  Mass  changes  at 
the  profile  and  farm  levels  are  most  sensitive  to 
variability  in  horizon  thickness  and  therefore  mass  of  soil 
which  influences  mass  of  organic  carbon.  This  is  especially 
the  case  for  Luvisols  because  mixing  of  the  thin  LFH  with 
subsurface  horizons  markedly  influences  mass  of  soil  in  the 
Ap  horizon. 

The  first  method  was  used  for  calculating  the  results 
presented.  It  seems  to  be  the  method  frequently  used  in  the 
literature.  In  many  cases,  however,  a  description  of  the 
method  used  is  missing. 


%  change= ( ( B-A)/A) 1 00 ;  where  A  is  original  carbon  content 
and  B  present  carbon  content;  negative  sign  means  loss, 
positive  sign  means  gain 
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4.3.2  Concentration  of  Organic  Carbon  in  Cultivated  and 
Native  A  Horizons 

Concentration  of  organic  carbon  in  native  surface 
horizons  was  highest  in  the  Gray  Luvisolic  zone(9.3%)  which 
contains  some  LFH  horizons  (Fig.  7).  Native  A  horizons  of 
the  Peace  River  area  include  some  Luvisols  accounting  for 
the  high  concentration  of  organic  carbon ( 8 . 3% ) ;  eliminating 
the  Luvisols  in  this  zone  lowers  the  average  concentration 
below  that  for  the  Black  soil  zone.  In  the  Chernozemic  Great 
Groups  there  was  an  expected  decline  in  concentration  of 
organic  carbon  in  the  native  A  horizons  with  increasing 
moisture  deficit  from  north  to  south  with  7.0%  for  the 
Black,  4.7%  for  the  Dark  Brown  and  2.9%  for  the  Brown  soil 
zone.  Average  concentration  for  the  Dark  Gray  soils  was 
3.3%. 

Upon  cultivation  the  concentration  in  organic  carbon 
decreased  with  the  greatest  absolute  difference  in  soils 
originally  highest  in  organic  carbon.  In  cultivated  A 
horizons  the  highest  concentration  of  organic  carbon  was 
present  in  the  Black  Chernozemic  soils(4.2%),  followed  by 
the  Dark  Gray(3.3%),  Peace  River(3.0%),  Dark  Brown(2.7%), 
Gray  Luvi sol ic ( 2 . 3% )  and  the  Brown  soils(1.4%). 

Differences  in  concentration  of  organic  carbon  between 
cultivated  and  native  A  horizons  were  highly  significant 
except  for  the  Dark  Gray  soils  which  might  have  resulted 
from  the  low  number  of  farms  sampled  in  this  zone. 
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Fig. 7  Average  Concentration  of  Organic  Carbon  in 

Cultivated  and  Native  A  and/or  LFH  Horizons  of  Soils 
in  Alberta. 
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4.3.3  Mass  of  Organic  Carbon  in  Cultivated  and  Native  A 
Horizons  and  Profiles 

Mass  of  organic  carbon  in  native  A  horizons  was  highest 
in  the  Black  soils  followed  by  the  DG,  PR,  GL ,  DB  and  BR 
soils  (Fig. 8);  in  the  native  profiles,  mass  follows  the  same 
distribution  except  that  PR  and  DG  soils  exchange  position. 
The  mass  of  OC  in  the  native  A  horizon  constitutes  47%,  58%, 
61%,  69%,  63%  and  55%  of  the  total  mass  of  organic  carbon  in 
the  profile  of  the  BR,  DB ,  BL,  DG,  GL  and  PR  soils. 

Mass  in  cultivated  A  horizons  and  profiles  was  lower 
than  in  native  samples.  Apart  from  the  DB  and  DG  soils,  mass 
of  organic  carbon  in  the  A  horizon  constitutes  a  higher 
proportion  of  the  total  mass  in  the  profile  than  in  the 
native  samples.  This  is  most  notable  for  the  Gray  Luvisols. 
Cultivated  A  horizons  contain  60%,  56%,  65%,  64%,  76%  and 
63%  of  the  total  mass  of  organic  carbon  in  the  profile. 

Differences  in  mass  of  organic  carbon  between 
cultivated  and  native  samples  were  highly  significant  for  A 
horizons  of  the  BR  and  DB  soils,  and  for  the  profiles  of  the 
BR,  DB  and  BL  soils;  differences  in  mass  were  significant 
for  the  DG  and  GL  profiles.  None  of  the  differences  for  A 
horizons  or  profiles  were  statistically  significant  for  the 
Peace  River  soils. 

At  the  individual  farm  level  some  positive  changes 
occurred:  21%, 13%,  33%, 43%  and  42%  of  the  sites  in  the  BR, 
DB,  BL,  GL  and  PR  sampling  zone  had  a  gain  in  mass  of 
organic  carbon  in  the  A  horizon.  Continuous  cropping  was 
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Sampling  Zone 


Ap  ^  Cultivated  profile  to  bottom 

]  Ah/LFH  of  B  horizon 

□  Native  Profile  to  bottom  of 
B  horizon 

Fig. 8  Average  Mass  of  Organic  Carbon  in  Cultivated  and 

Native  A  and/or  LFH  Horizons  and  Profiles  of  Soils 
in  Alberta. 


53 


practiced  on  34% , 1 00% , 60% , 34%  and  60%  of  these  sites  in  the 
above  order  of  the  sampling  zones.  Lay  farming  was  practiced 
on  20%  in  the  Black  and  33%  in  the  Gray  Luvisolic  zone.  The 
remainder  is  constituted  by  fallow  rotations.. 

A  mass  gain  in  the  profile  occurred  on  14%, 13%, 14%  and 
33%  of  the  sites  in  the  BR,  BL ,  GL  and  PR  zones. 
Continuously  cropped  sites  constituted  50%  of  these  sites  in 
the  BR,  BL  and  PR  zones;  50%  of  these  sites  in  the  GL  zone 
were  under  lay  farming;  the  remaining  percentages  were  under 
fallow  rotations. 

4.3.4  Changes  in  Concentration  and  Hass  of  Organic  Carbon  in 
A  Horizons 

The  change  in  mass  was  less  than  in  concentration 
(Fig. 9)  except  for  the  DG  soils  where  it  was  highly 
significantly  greater.  Differences  were  highly  significant 
for  the  BR,  DB  and  GL  soils;  they  were  significant  for  the 
BL  and  PR  soils.  Numerical  values  are  listed  in  Appendix 
II  .  1  to  II .6. 

4.3.5  Average  Thickness  of  Cultivated  and  Native  A  Horizons 

Apart  from  the  Dark  Gray  Chernozemic  soils,  cultivated 
A  horizons  were  thicker  than  native  A  horizons  (Fig. 10) 
suggesting  mixing  of  subsurface  horizon  material  into  the  A 
horizon.  This  trend  was  especially  apparent  for  the  Gray 
Luvisolic  and  Peace  River  soils.  Due  to  great  variability  in 
horizon  depths,  however,  none  of  the  differences  in  horizon 
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Fig. 9  Average  Loss  in  Concentration  and  Mass  of  Organic 
Carbon  from  A  and/or  LFH  Horizons  of  Soils  in 
Alberta  following  Cultivation. 
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Fig.  10  Average  Thickness  of  Cultivated  and  Native 
LFH  Horizons  of  Soils  in  Alberta. 
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pairs  was  statistically  significant.  Similar  variability  in 
horizon  depths  was  reported  by  Thiessen  et  al.(1982). 
Caldwell  et  al.(1939)  found  considerable  variability  in 
horizon  depths  of  duplicate  profiles. 

4.3.6  Average  Apparent  Bulk  Density  of  Cultivated  and  Native 
A  Horizons 

Apparent  bulk  density  of  A  horizons  increased  through 
cultivation  (Fig. 11).  Differences  for  the  Gray  Luvisolic  and 
the  Peace  River  soils  were  highly  significant;  differences 
were  significant  for  the  Brown  and  Black  zones.  There  was  no 
statistically  significant  difference  for  the  Dark  Brown  and 
Dark  Gray  sampling  zones. 

4.3.7  Average  pH  of  Cultivated  and  Native  A  Horizons 

A  trend  toward  higher  pH  which  might  have  been  due  to 
horizon  mixing  was  observed  in  5  zones  (Fig . 1 2 ) . Thi s  was 
only  significant  for  the  Brown  and  Gray  Luvisolic  Great 
Groups . 

4.3.8  Distribution  of  Organic  Carbon,  Apparent  Bulk  Density 
and  pH  with  Depth  in  Cultivated  and  Native  Profiles 

The  distribution  of  apparent  bulk  density,  pH  and 
concentration  of  organic  carbon  with  depth  (Fig.s  13  to  18) 
shows  that  changes  in  these  properties  upon  cultivation  also 
occurred  in  subsurface  horizons.  Profile  depths  to  the 
bottom  of  the  B  horizon  were  similar  in  cultivated  and 
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Fig.  1  1 


Average  Apparent  Bulk  Density  of  Cultivated  and 
Native  A  and/or  LFH  Horizons  of  Soils  in  Alberta. 
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Fig. 12  Average  pH  in  Water  of  Cultivated  and  Native  A 
and/or  LFH  Horizons  of  Soils  in  Alberta. 
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Fig. 13  Distribution  of  Organic  Carbon,  Apparent  Bulk 

Density  and  pH  with  Depth  in  Cultivated  and  Native 
Soils  in  the  Brown  Sampling  Zone. 
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Fig. 14  Distribution  of  Organic  Carbon,  Apparent  Bulk 

Density  and  pH  with  Depth  in  Cultivated  and  Native 
Soils  in  the  Dark  Brown  Sampling  Zone. 
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Fig. 15  Distribution  of  Organic  Carbon,  Apparent  Bulk 

Density  and  pH  with  Depth  in  Cultivated  and  Native 
Soils  in  the  Black  Sampling  Zone. 
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Fig. 16  Distribution  of  Organic  Carbon,  Apparent  Bulk 

Density  and  pH  with  Depth  in  Cultivated  and  Native 
Soils  in  the  Dark  Gray  Sampling  Zone. 
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Fig.  17  Distribution  of  Organic  Carbon,  Apparent  Bulk 

Density  and  pH  with  Depth  in  Cultivated  and  Native 
Soils  in  the  Gray  Luvisolic  Sampling  Zone. 
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Fig. 18  Distribution  of  Organic  Carbon,  Apparent  Bulk 

Density  and  pH  with  Depth  in  Cultivated  and  Native 
Soils  in  the  Peace  River  Sampling  Zone. 
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native  profiles  except  for  the  Gray  Luvisolic  soils  where 
the  cultivated  profile  averaged  10  cm  deeper  than  the  native 
profile.  In  all  sampling  zones  concentration  of  organic 
carbon  decreased  with  depth,  apparent  bulk  density 
increasing  simultaneously.  With  the  exception  of  the  Gray 
Luvisolic  profiles,  pH  increased  with  depth.  It  was  higher 
in  the  cultivated  than  in  the  native  profiles  with  the 
exception  of  the  Dark  Gray  soils.  Organic  carbon  content  in 
the  cultivated  B  horizon  was  lower  than  in  the  native  B 
horizon  in  the  Brown,  Black  and  Gray  Luvisolic  soils.  It  was 
higher  in  the  cultivated  than  in  the  native  B  horizon  for 
the  Dark  Brown,  Dark  Gray  and  Peace  River  soils.  Changes  in 
apparent  bulk  density  were  opposite  to  changes  in 
concentration  of  organic  carbon.  Differences  in 
concentration  of  organic  carbon  in  the  B  horizon  of  the  Dark 
Brown  soils  and  pH  in  the  B  horizon  of  the  Brown  soils  were 
highly  significant.  None  of  the  remaining  differences  in 
properties  between  cultivated  and  native  subsurface  horizons 
were  statistically  significant. 


4.3.9  Moisture  Retention  in  Cultivated  and  Native  A  Horizons 

Moisture  retention  at  -33  and  -1500kPa  in  cultivated  A 
horizons  was  lower  than  in  native  A  horizons  in  all  sampling 
zones  (Fig. 19).  Differences  were  highly  significant  for  the 
Dark  Brown,  Black,  and  Peace  River  zones, 
statistically  significant  in  the  remaining  zones. 
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Fig. 19  Average  Moisture  Retention  at  -33  and  -1500 
Cultivated  and  Native  A  and/or  LFH  Horizons 
in  Alberta. 
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In  all  cases,  moisture  retention  decreased  more  at  the 
lower  potential  than  at  the  higher  potential.  The  result  was 
that  the  difference  between  the  two,  used  as  an  indicator 
for  %  available  moisture  capacity,  was  higher  in  cultivated 
than  in  native  A  horizons.  Differences  in  %  available 
moisture  were  highly  significant  for  the  Brown,  Dark  Brown 
and  Black  soils,  but  not  statistically  significant  for  the 
remaining  zones. 

Distribution  of  moisture  retention  at  -33  and  -1500kPa 
over  the  sampling  zones  followed  the  same  pattern  as  the 
distribution  of  the  mass  of  organic  carbon  in  Figure  8. 
Moisture  retention  at  -33  and  -1500kPa  was  positively 

f 

correlated  with  concentration  of  organic  carbon  and  texture 
of  the  parent  material;  it  was  negatively  correlated  with 
apparent  bulk  density  (Tables  6  and  7).  Available  moisture 
capacity  in  native  and  cultivated  A  horizons  on  the  other 
hand  increased  from  south  to  north  in  the  sampling  zones 
(Fig. 20)  . 

4.4  Discussion 

Significant  losses  in  organic  carbon  on  a  concentration 
and  mass  basis  have  occurred  in  the  six  sampling  zones  (Fig. 
7,8,9).  The  data  show  that  changes  not  only  affected  surface 
horizons  but  the  complete  profile  (Fig.  8,13  to  18). 

Similar  trends  were  reported  by  Newton  et  al. ( 1945) , 
Dormaar  (1979),  Dormaar  and  Pittman  (1980),  McGill  et 
al. ( 1981 ) ,  Voroney  et  al.(1981)  and  Campbell 


and  Souster 


Table  6  Significant  ( p<0 . 05 ) Corre 1  a 1 1  on  Coefficients  Between  Moisture  Retention  and  Management  and  Soil  Variables  In 

Cultivated  A  Horizons 
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Table  7  Significant  ( p<0 . 05 ) Cor re  1  at  1  on  Coefficients  Between  Moisture  Retention  and  Soil  Variables  in  Native 

Hor i zons 
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Fig. 20  Average  Available  Moisture  Capacity  in  Cultivated 
and  Native  A  and/or  LFH  Horizons  of  Soils  in 
Alberta . 
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(1982)  but  with  different  magnitude  of  the  changes.  Losses 
in  concentration  of  organic  carbon  for  the  Brown  soils  are 
higher  than  those  reported  by  Dormaar ( 1 97 9 )  and  Campbell  and 
Souster (  1 982 )  (see  Table  1).  Losses  for  the  Dark  Brown  soils 
are  lower  than  comparable  values  for  fallow  rotations 
reported  by  the  same  authors,  but  somewhat  higher  than 
losses  for  a  continuous  crop  rotation  reported  by  Dormaar 
and  Pi ttman ( 1 980 ) .  Different  magnitudes  of  the  changes  might 
be  due  to  the  different  numbers  of  sites  on  which  various 
results  are  based;  as  well  as  different  crop  rotations  and 
years  of  cultivation.  Results  of  Dormaar ( 1 979 ) ,  Dormaar  and 
Pi ttman ( 1 980 )  and  Campbell  and  Souster ( 1 982 )  were  derived 
using  data  from  one,  one,  and  three  sites,  respectively, 
with  mainly  crop-summer  fallow  rotations.  Results  here  are 
based  on  14  sites  in  the  Brown  and  16  sites  in  the  Dark 
Brown  zone  integrating  fallow  rotations  and  continuous 
cropping  with  various  cultivation  times. 

Comparing  literature  results  to  data  for  farms  from  the 
present  study  with  similar  rotation  and  cultivation  time, 
however,  yields  closer  agreement.  For  example,  two  of  the  14 
sites  in  the  Brown  sampling  zone  and  three  of  the  16  in  the 
Dark  Brown  zone  had  losses  that  approximated  those  reported 
in  the  literature.  Losses  reported  by  Campbell  and 
Souster ( 1 982 )  for  the  Dark  Brown  soils  lie  within  the  range 
of  losses  for  11  fallow  rotations  of  the  16  Dark  Brown 


sites. 
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Mass  losses  in  the  A  horizons  are  only  about  1/3  to  2/3 
of  those  estimated  by  McGill  et  al.(1981),  and  1/3  the 
average  loss  from  3  slope  positions  at  a  Black  Chernozemic 
site  from  Saskatchewan  (Voroney  et  al.,1981;  Table  2).  The 
differences  are  again  probably  due  to  the  present  larger 
data  base,  and  perhaps  also  to  crop  rotation.  None  of  the 
present  results  for  A  horizons  on  an  individual  farm  basis 
(Appendix  II. 3)  are  similar  to  those  of  Voroney  et  al.  (  1981 ) 
(Table  2 ) . 

McGill  et  al. ( 1981  )  predicted  the  lowest  proportional 
loss  for  the  Gray  Luvisols  despite  their  faster 
decomposition  rate  and  lower  initial  amount  of  organic 
matter.  The  present  results  suggest  that  the  Black  soils 
with  the  slowest  decomposition  rate  and  highest  initial 
amount  had  the  lowest  proportional  loss.  Estimates  by  McGill 
et  al. (  1981 )  of  initial  amounts  are  also  consistently  lower 
than  present  results.  This  creates  greater  proportional 
losses  because  of  the  smaller  denominator.  Their  data  used 
as  the  original  amount  of  organic  carbon  were  averages  from 
CanSIS  which  included  some  Ap  as  well  as  Ah  horizons  and 
used  literature  average  densities  for  conversion  to  mass. 
Their  calculations  of  mass  changes  did  not  adequately  treat 
concomitant  density  increase. 

On  the  other  hand,  there  is  agreement  with  data  of 
Newton  et  al.(1945)  for  the  Black  and  Gray  Luvisolic 
sampling  zones  although  average  cultivation  time  has  nearly 
tripled  and  Newton’s  samples  were  taken  from  crop-fallow 
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rotations.  Newton  et  al.(1945)  reported  mass  changes  of  -18% 
and  -27%  for  the  Black  and  Gray  Luvisolic  zones  in  1945 
compared  to  observed  changes  of  -17%  and  -27%  based  on 
sampling  in  1982;  losses  from  the  Brown  and  Dark  Brown  A 
horizons  are  higher  than  in  Newton's  reports.  The  number  of 
sampling  locations  was  similar  in  both  studies;  depth 
increments  of  0-15  cm  used  by  Newton  et  al.(1945)  are 
similar  to  the  present  average  horizon  depths  of  cultivated 
A  horizons;  different  methodologies  for  carbon  determination 
(dry  combustion  vs.  wet  oxidation)  were  used,  however;  and 
bulk  densities  used  in  calculating  mass  were  not  reported  by 
Newton  et  al.(1945). 

Two  possible  patterns  of  mass  changes  in  organic  carbon 
may  have  occurred  in  the  A  horizons:  i)losses  in  the  Black 
and  Gray  Luvisolic  zones  approached  a  steady-state  by  1945 
and  have  not  changed  much  since,  whereas  losses  continued  in 
the  Brown  and  Dark  Brown  zones;  or  ii)organic  carbon 
continued  to  decrease  after  Newton's  investigations  but,  at 
some  point  in  time,  increased  again,  due  perhaps  to 
introduction  of  different  management  practices  on  a  more 
extensive  basis  in  the  BL  and  GL  zones  than  in  the  DB  and 
BR.  The  overall  change  in  organic  carbon  could  therefore  be 
similar  now  to  that  reported  40  years  ago  for  the  BL  and  GL 
soils.  Organic  matter  levels  for  the  BR  and  DB  soils  are 
still  below  those  reported  by  Newton  at  al.(1945).  The 
second  scenario  suggests  higher  organic  carbon  losses  would 
have  been  found  during  some  time  between  1945  and  1982  than 
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were  observed  in  1982. 

Mass  losses  from  the  profiles  of  the  Black  soils  are 
half  of  the  loss  (21%  vs.  42%)  calculated  from  data  of 
Voroney  et  al .( 1 98 1 )( Table  2),  but  they  had  only  one  sample 
site.  Their  reported  losses  are  within  the  range  of  losses 
for  crop-crop-fallow  and  crop-fallow  rotations  for  4  farms 
of  the  Black  zone  in  the  present  study  (Appendix  II. 3). 

Mass  losses  in  the  subsurface  horizons  are  up  to  6 
times  higher  than  those  reported  by  Newton  et  al.(1945)  for 
the  15-30  cm  depth.  Despite  the  difference  in  depth 
increments  and  absence  of  bulk  density  data  in  reports  by 
Newton  et  al.(1945),  there  is  a  suggestion  that  organic 
carbon  in  subsurface  horizons  has  continued  to  decrease, 
even  if  changes  in  surface  horizons  might  have  stopped  or 
altered  direction. 

The  two  mechanisms  causing  organic  carbon  losses  are 
decomposition  and  erosion.  Erosion  was  not  measured.  No 
reduction  in  profile  depths  seemed  to  have  occurred  (Fig. 13 
to  18).  One  has  to  conclude  that  no  measurable  net  amount  of 
soil  was  removed  by  erosion.  Assessment  of  the  relative 
magnitudes  of  the  two  processes  is  not  possible.  The 
variability  in  profile  thicknesses  causes  such  a  measure  of 
erosion  to  have  a  low  resolution. 


Lower  losses  in  mass  of  organic  carbon  as  opposed  to 
concentration  of  organic  carbon  (Fig. 9)  may  be  due  to  two 
factors.  First  is  mixing  of  subsurface  horizon  material  into 
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the  A  horizon  during  cultivation.  This  dilutes  the 
concentration  of  organic  carbon.  A  more  negative  change  in 
concentration  of  organic  carbon  results,  but  with  no 
corresponding  change  in  mass.  Changes  in  A  horizon 
thicknesses  (Fig. 10)  were  greatest  in  soils  of  the  Gray 
Luvisolic  and  Peace  River  zones  which  also  have  the  highest 
losses  on  a  concentration  basis.  Further,  'thickness  of 
cultivated  B  horizons  tended  to  decrease  (Fig. 13  to  18) 
which  is  consistent  with  its  incorporation  into  Ap  horizons. 
Second  is  an  inverse  relationship  between  changes  in 
concentration  of  organic  carbon  and  apparent  bulk  density 
(Fig. 21).  Russell ( 1 960 )  and  Davidson  et  al.(1967)  reported 
an  inverse  relationship  between  concentration  of  organic 
carbon  and  bulk  density.  This  is  also  evident  from  Fig. 11 
and  Figures  13  to  18.  Correlation  coefficients  between 
apparent  bulk  density  and  concentration  of  organic  carbon  in 
cultivated  A  horizons  were  -.5  for  the  Brown,  Dark  Brown  and 
Black  soils  and  -.6  for  the  Gray  Luvisols  in  the  present 
study.  These  coefficients  were  not  statistically 
s ign i f icant . 

Deepening  of  A  horizons,  reduction  in  concentration  of 
organic  carbon,  and  increase  in  apparent  bulk  density  are 
linked.  The  mass  of  soil  in  cultivated  A  horizons  is  greater 
than  in  native  A  horizons.  The  increase  in  soil  mass  exceeds 
the  increase  in  depth.  Therefore,  a  larger  amount  of  B 
horizon  material  entered  the  A  horizon  than  would  be 
suggested  by  changes  in  horizon  thickness.  The  associated 


Change  in  Apparent  Bulk  Density  (%) 
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Change  in  Concentration  of  Organic  Carbon  (%) 


Fig. 21  Relationship  Between  Change  in  Apparent  Bulk  Density 
and  Change  in  Concentration  of  Organic  Carbon  in  A 
Horizons  of  Soils  in  Alberta  (r=0.47;  p<0.0001; 
n  =  64 )  . 
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increase  in  bulk  density  has  prevented  the  extent  of 
dilution  from  being  obvious. 

Concern  about  changes  in  concentrat ion  or  changes  in 
mass  of  organic  carbon  will  depend  on  whether  qualitative 
(e.g.  aggregation)  or  quantitative  (e.g.  nutrient  supply) 
aspects  of  organic  matter  are  of  interest.  With  regard  to 
maintenance  of  productivity  both  are  equally  important. 


Available  moisture  capacity  in  cultivated  A  horizons 
was  higher  than  in  native  A  horizons  in  all  sampling  zones 
because  moisture  retention  decreased  more  at  -1500kPa  than 
at  -33kPa.  Moisture  retention  at  -33kPa  and  -1500kPa 
increased  with  increasing  mass  of  organic  carbon.  This 
suggests  that  organic  matter  influences  moisture  retention 
both  at  low  and  high  potentials,  but  apparently  more  at  the 
lower  than  at  the  higher  potential.  Such  a  conclusion  was 
also  reached  by  De  Jong  et  al.(1983). 

Clay  and  bulk  density  are  also  positively  correlated 
with  moisture  retention  (Hutcheon , 1 942 ;  Hill  and  Sumner, 
1967;  Shaykewich  and  Zwar ich , 1 968 ;  Reeve  et  al.,1973;  De 
Jong  et  al.,1983).  Shaykewich ( 1 980 )  concluded  that  particle 
size  distribution  was  more  important  in  its  influence  on 
available  moisture  capacity  than  organic  matter.  The 
tendendy  for  clay  contents  to  increase  from  south  to  north 
throughout  the  province  may  be  related  to  the  increase  in 
available  moisture  capacity  (Fig. 20).  No  particle  size 


determination  was  conducted  for  the  individual  samples; 
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texture  of  the  parent  material  was  estimated  by  hand. 

Actual  moisture  available  under  field  conditions,  is 
additionally  influenced  by  pore  size  distribution  and 
infiltration  rate.  Both  of  these  properties  depend  on 
texture-organic  matter-bulk  density  relationships.  The 
linking  factor  is  aggregation.  Organic  matter  improves 
aggregation  and  root  penetration;  it  thus  allows  a  larger 
volume  of  soil  to  be  explored  by  roots  and  potentially  more 
moisture  to  be  taken  up  by  plants  ( Allison , 1 973 ) . 

Determination  of  available  moisture  capacity  on 
disturbed  samples  will  therefore  not  completely  describe  the 
effect  of  changes  in  organic  carbon  on  moisture  regimes  in 
the  field. 


4.5  Conclusions 

Significant  losses  in  concentration  and  mass  of  organic 
carbon  in  A  horizons  and  profiles  occurred  through 
cultivation.  Losses  in  concentration  of  organic  carbon  were 
significantly  higher  than  losses  in  mass  of  organic  carbon 
in  most  cases.  This  has  been  attributed  to  the  combined 
effects  of;  i)  dilution  due  to  mixing  of  subsurface  horizon 
material  into  the  original  A  horizon;  and  ii)  an  inverse 
relationship  between  changes  in  apparent  bulk  density  and 
changes  in  concentration  of  organic  carbon.  Both  of  these 
mechanisms  are  linked. 

Depending  on  the  sampling  zone,  variable  gains  in  mass 
of  organic  carbon  occurred  on  13%  to  43%  of  the  farms  for  A 
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horizons  and  on  13%  to  33%  of  the  farms  for  the  profile. 

Thickness  of  A  horizons  increased  due  to  cultivation 
except  in  the  Dark  Gray  soils.  None  of  the  differences  in 
horizon  pairs  was  statistically  significant  which  was 
attributed  to  great  variability  in  horizon  thickness. 

Apparent  bulk  density  in  A  horizons  increased  due  to 
cultivation  in  all  sampling  zones.  The  increase  was  not 
statistically  significant  in  the  Dark  Brown  and  Dark  Gray 
sampling  zones. 

No  reduction  in  profile  depths  occurred  due  to 
cultivation.  Losses  of  organic  carbon  therefore  were 
attributed  mainly  to  decomposition.  No  measurement  of 
erosion  was  conducted  to  confirm  this  conclusion. 

Organic  carbon  changes  influenced  moisture  retention 
more  at  lower  than  at  higher  potentials.  Consequently 
available  moisture  capacity  increased  in  cultivated  A 
horizons . 

The  magnitude  of  the  change  in  organic  carbon  depends 
on  the  method  of  calculation.  The  method  used  in  this  study 
corresponds  to  that  used  in  the  literature  and  normally 
yields  the  most  negative  changes. 

Losses  reported  in  the  literature  correspond  to  two  of 
the  14  sites  in  the  Brown  sampling  zone  and  three  of  the  16 
sites  in  the  Dark  Brown  sampling  zone,  if  crop  rotations  and 
cultivation  time  were  similar.  Although  the  literature 
values  correspond  to  some  individual  farms,  they  do  not 
appropriately  represent  the  average  for  any  soil  zone.  A 
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larger  number  of  sites  is  needed  than  has  been  used  in  the 
past.  Comparison  with  data  from  Newton  et  al.(1945)  suggests 
that  losses  may  have  ceased  in  the  Black  and  Gray  Luvisolic 
sampling  zones. 


_  J 


5.  Regional  and  Management  Variables  Affecting  Changes  in 

Soil  Organic  Carbon 

5 . 1  Introduction 

The  objective  of  this  section  is  to  examine  the  effects 
of  regional  and  management  variables  on  changes  in  organic 
carbon  presented  in  the  previous  chapter.  Management  data 
used  cover  the  period  from  1982  to  1967.  Many  missing  data 
between  1967  and  1907  made  regression  analyses  for  that 
period  impossible.2 


5.2  Methods  of  Statistical  Analysis 

The  overall  dependence  of  a  criterion  variable  such  as 
present  concentration  of  organic  carbon  on  a  single  or  set 
of  independent  (predictor)  variables  (management  data,  soil 
properties)  was  assessed  by  simple  linear  and  stepwise 
forward  multiple  linear  regression  techniques.  Statistical 
analysis  was  performed  using  SPSSx  subprograms  Scattergram 
and  Regression.  Only  farms  with  complete  data  sets  were 
analyzed  statistically.  The  Dark  Gray  sampling  zone, 
represented  by  only  two  farms,  was  eliminated. 


2 

Data  codes  are  listed  in  Appendix  I.  Ascending  numerical 
values  for  crop  rotation  represent  lower  fallow  frequency 
and  more  continuous  cropping  as  numbers  increase;  ascending 
numbers  for  organic  residue  represent  quantitatively  (e.g. 
straw  chopped)  and  qualitatively  (e.g. manure)  higher  return. 
Tillage  intensity  and  input  of  N  and  P  increase  with 
ascending  numbers. 
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5.3  Results 

5.3.1  Characteristics  of  Management  Practices  throughout 
Alberta 

The  proportion  of  various  management  practices  in  the 
sampling  zones  is  summarized  in  Table  8.  Percentages  are  of 
the  total  number  of  farms  for  which  data  were  available. 

Fields  in  the  Brown  sampling  zone  had  the  lowest 
average  cultivation  time  of  38  years.  No  data  were  available 
for  three  farms  and  four  out  of  the  remaining  11  fields  had 
been  cultivated  for  only  three,  six,  ten  and  18  years, 
respectively.  Fields  in  the  Black  and  Dark  Brown  zone  had 
similar  cultivation  time  of  around  60  years;  data  for  one 
farm  were  missing  in  the  Black  zone.  The  Peace  River  area 
was  under  cultivation  longer  than  the  Gray  Luvisolic  zone; 
data  from  three  farms  were  missing  in  the  Peace  River  area. 

On  a  provincial  level,  summerf allowing  in  combination 
with  stubble  mulching  was  the  most  common  practice,  being 
applied  to  51%  of  the  fields  sampled.  Some  of  these  fields 
were  cultivated  with  a  noble  blade  to  reduce  tillage 
intensity.  Continuous  cropping  on  the  other  hand  had  a 
frequency  of  44%.  Mineral  N  and  P  were  applied  to  86%  of  the 
fields  during  some  time  in  the  15  years  previous  to 
sampling;  20%  of  the  total  fields  received  fertilizer  only 
during  five  years  prior  to  sampling. 

Within  the  sampling  zones,  fields  of  the  Brown  zone  had 
the  highest  proportion  of  summer  fallowing  and  stubble 


TABLE  8:  RELATIVE  FREQUENCIES  OF  MANAGEMENT  PRACTICES  ON  FARMS  SAMPLED  IN  ALBERTA  (1982-1967) 
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mulching  followed  by  the  Dark  Brown  zone;  in  both  zones  use 
of  the  noble  blade  was  not  uncommon.  More  continuous 
cropping  was  practiced  on  fields  in  the  Dark  Brown  zone  than 
in  the  Brown  zone. 

Whereas  nearly  70%  of  the  fields  sampled  in  the  Dark 
Brown  zone  had  regular  inputs  of  mineral  N  and  P  since  1967, 
only  18%  in  the  Brown  zone  did  so.  In  the  Brown  zone  46%  of 
the  fields  had  not  received  mineral  fertilizer  until  five 
years  prior  to  sampling  time.  The  highest  relative  frequency 
of  continuous  cropping,  manuring,  amount  of  straw  returned, 
N  and  P  applications  and  conventional  tillage  was  on  fields 
of  the  Black  sampling  zone.  Continuous  cropping  and 
application  of  mineral  fertilizer  had  second  highest 
frequency  in  the  Peace  River  area  after  the  Black  sampling 
zone.  Continuous  cropping  and  lay  farming  were  not  uncommon 
in  the  Gray  Luvisolic  zone  and  mineral  N  and  P  were  used 
more  frequently  than  in  the  Brown  zone.  Green  manuring  had 
frequencies  of  14%  and  11%  in  the  Gray  Luvisolic  and  Peace 
River  zones,  respectively. 

Fertilizer  N  input  started  some  time  between  1947  and 
1967  (Table  9).  Annual  inputs  to  fallow  rotations  were  very 
low.  Mineral  P  was  not  applied  until  some  time  after  1967 
(Table  10).  Simultaneously,  continuous  cropping  became  more 
frequent  after  1967  entailing  a  considerable  increase  in 
mineral  N  and  P  input  (Table  9  and  10).  The  highest  average 
annual  input  of  N  between  1982  and  1977  was  in  the  Black  and 
Peace  River  zones  under  continuous  cropping.  The  highest 


TABLE  9  AVERAGE  ANNUAL  N  INPUTS  (kg-ha1)  TO  CEREAL  CROPS  UNDER  VARIOUS  ROTATIONS  BETWEEN  1982  AND  1907. 
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average  annual  P  input  during  the  same  period  was  to 
crop-crop-summer  fallow  rotations  in  the  Gray  Luvisolic  zone. 
In  the  remaining  areas,  however,  N  and  P  input  were  many 
times  higher  under  continuous  cropping  than  under  fallow 
rotations. 

Cereal  yields  seem  to  have  increased  concomitantly  with 
fertilizer  input  (Table  11).  This  is  consistent  with'yield 
trends  for  the  three  prairie  provinces  reported  by  Ridley 
and  Hedlin  (1980).  Highest  average  annual  grain  yields  were 
achieved  under  continuous  cropping;  under  this  management 
the  Black  sampling  zone  had  the  highest  average  yield 
amongst  the  sampling  zones. 


5.3.2  Relationship  between  Crop  Rotation  and  Change  in 
Organic  Carbon 

Losses  in  concentration  and  mass  of  organic  carbon  were 
consistently  lower  under  continuous  cropping  than  under 
fallow  rotations  in  the  Brown,  Dark  Brown  and  Black  zones 
(Table  12).  Higher  losses  under  continuous  cropping  in  the 
Gray  Luvisolic  and  Peace  River  sampling  zones  were  probably 
due  to  the  presence  of  LFH  horizons  amongst  the  native 
horizons . 


5.3.3  Simple 

Simple 

the  dependen 


Linear  Regression  Models 
linear  regressions  were 
t  variables: 


performed  for  each  of 


1.  concentration  of  organic  carbon  in  cultivated  A  horizons; 


TABLE  11  AVERAGE  ANNUAL  YIELDS  (Y)  (Mg.ha'‘)OF  CEREAL  CROPS  UNDER  VARIOUS  ROTATIONS  BETWEEN  1982  AND  1907. 
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TABLE  12  AVERAGE  CHANGE  IN  CONCENTRATION  AND  MASS  OF  ORGANIC  CARBON(OC)  IN  A  HORIZONS  UNDER  SUMMERFALLOW  ( SF ) AND  UNDER 

CONTINUOUS  CROPPING  (CC) 
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2.  change  in  concent rat i on  of  organic  carbon  in  the  A 
horizon ; 

3.  change  in  mass  of  organic  carbon  in  the  A  horizon; 
with  the  predictor  variable  years  of  cultivation. 

The  purpose  was  to  examine  whether  there  was  a  defined 
relationship  between  any  of  the  independent  variables  and 
the  criterion  variable;  and  to  determine  the  direction  of 
this  relationship.  The  overall  pattern  of  the  scatter  was 
used  as  an  indicator  for  any  systematic  relationship  among 
the  data  represented. 


Relationship  between  Change  in  Organic  Carbon  and  Time 

On  a  provincial  level,  a  weak  but  positive  correlation 
existed  between  present  concentration  of  organic  carbon  and 
time  (Appendix  1 1 1 .  1  )  . 

On  a  sampling  zone  level  the  only  statistically 
significant  trend  was  in  the  Brown  sampling  zone,  where  31% 
(p<0.05)  of  the  variability  in  change  in  concentration  of 
organic  carbon  was  explained  by  years  of  cultivation,  the 
direction  being  negative  (Appendix  III. 2).  Data  points 
showed  a  curvilinear  trend  (Fig. 22). 

The  change  in  organic  carbon  with  time  has  been 
described  as  a  negative  curvilinear  relationship  (Henin  and 
Dupuis, 1945;  Woodruf f , 1 949 ) .  Initially,  organic  carbon 
decreases  rapidly  with  the  slope  levelling  off  afterwards. 
Apart  from  the  Brown  zone  (Fig. 22)  none  of  the  scattergram 
patterns  indicated  such  a  relationship.  The  Black  sampling 
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Change  in  Concentration  of 
Organic  Carbon i(%) 
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Brown  Sampling  Zone 


Years  of  Cultivation 


Fig. 22  Relationship  Between  Change  in  Concentration  of 

Organic  Carbon  and  Time  in  the  Brown  Sampling  Zone 
(r=-0.56;  p=0.04;  n=11). 


zone 
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(Fig. 23)  was  chosen  to  represent  trends  in  the 
remaining  sampling  zones.  There  did  not  seem  to  be  any 
systematic  relationship  with  time  among  the  data 
represented . 

In  the  Gray  Luvisolic  zone  a  positive  correlation  of 
+.58  existed  between  change  in  mass  of  organic  carbon  and 
years  of  cultivation  at  a  significance  level  of  8%  (Appendix 
III  .3)  . 

Influence  of  Management  Practices  on  Change  in  Organic 

Carbon  With  Time 

A  second  set  of  simple  linear  regressions  was  performed 
for  the  dependent  variables  described  and  cultivation  time 
as  a  predictor  variable  when  different  kinds  of  crop 
rotations  were  practiced  and  different  kinds  of  organic 
residues  returned.  Crop  rotations  and  organic  residues  were 
as  described  in  Table  8  and  as  coded  in  Appendix  I  (File 
Cropinfo).  The  objective  was  to  determine  if  prediction 
accuracy  of  change  in  organic  carbon  with  time  could  be 
improved . 

On  a  province  level,  13%  of  the  variability  (n=26)in 
present  concentration  and  change  in  concentration  of  organic 
carbon  was  explained  by  years  of  cultivation  if  continuous 
cropping  was  practiced.  The  slope  was  positive  (Appendix 
III. 4).  When  straw  was  removed  which  was  the  case  for  6 
farms,  80%  of  the  variability  in  present  concentration  and 
change  in  concentration  of  organic  carbon  was  explained  by 


' 


Change  in  Concentration  of 
Organic  Carbon  (%) 
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Black  Sampling  Zone 


Fig. 23  Relationship  Between  Change  in  Concentration  of 

Organic  Carbon  and  Time  in  the  Black  Sampling  Zone 
(  r  =  -0 . 05 ;  p=0.43;  n=14). 
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years  of  cultivation.  Continuous  cropping  was  simultaneously 
practiced  on  50%  of  these  farms,  however. 

On  a  sampling  zone  level,  four  statistically 

significant  trends  were  obvious: 

1.  Brown  zone:  In  stubble  mulch  systems  41%  of  the 
variability  in  the  change  in  mass  of  organic  carbon  was 
explained'  by  cultivation  time  (Appendix  III. 5).  Stubble 
mulching  in  combination  with  summer fallowi ng  was  common. 
Although  not  statistically  significant,  a  negative 


relationship  existed 

between 

change  in  concentration 

of 

organic 

carbon  and 

cultivation  time  if  summerf allowing 

was 

practiced  (Fig. 24). 

2 . Dark 

Brown  zone: 

If  straw  ’ 

was  removed  (three  farms), 

99% 

of  the 

variability 

(positive 

slope)  in  mass  change 

of 

organic 

carbon  was 

explained 

by  cultivation  time  (Appen 

dix 

III. 6).  Two  of  the  three  farms,  however,  had  practiced 
continuous  cropping  for  15  years.  The  third  farm  site  had 
been  under  a  crop-crop-fallow  rotation  for  60  years 
(1917-1977),  but  the  field  was  put  under  pasture  in  the  5 
years  previous  to  sampling  which  seems  to  have  had  an  effect 
on  rebuilding  organic  matter.  The  high  correlation  obtained 
has  to  be  interpreted  with  great  care,  however,  because  only 
three  data  points  were  present.  This  reduces  the  degrees  of 
freedom  to  a  minimum  and  can  result  in  apparently  high 
correlation  of  unrelated  populations. 

3.  Gray  Luvisolic  zone:In  crop-crop-fallow  rotations  97%  of 
the  variability  in  present  concentration  of  organic  carbon 


. 


Change  in  Concentration  of 
Organic  Carbon  (o/0) 
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Brown  Sampling  Zone 
Summer  Fallow 


Years  of  Cultivation 


Fig. 24  The  Influence  of  Summer  fallowing  on  the  Relationship 
Between  Change  in  Concentration  of  Organic  Carbon 
and  Time  in  the  Brown  Sampling  Zone  (r=-0.56; 
p=0 . 07 ;  n  =  8 ) . 
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was  explained  by  cultivation  time.  The  slope  was  slightly 
negative  (Appendix  III. 8).  This  was  the  case  for  only  three 
farms  and  the  high  correlation  is  of  little  meaning. 

4. Peace  River  zone:  When  summer  fallowing  was  practiced  95% 
of  the  variability  in  change  in  concentration  of  organic 
carbon  was  explained  by  cultivation  time.  The  slope  was 
positive  (Appendix  III. 9).  Again,  only  four  cases  were 
present  causing  problems  as  mentioned  for  the  Dark  Brown  and 
Gray  Luvisolic  zone. 

No  statistical  significance  existed  between  the  change 
in  concentration  of  organic  carbon  and  time  in  the  Black 
zone  when  continuous  cropping  was  practiced  (Fig. 25). 

Significance  and  correlation  between  change  in 
concentrat ion  of  organic  carbon  and  cultivation  time  were 
not  improved  by  taking  management  practices  into  account 
(statistics  at  bottom  of  Fig. 24  and  25).  The  direction  of 
the  relationship,  however,  remained  consistent.  Yet,  much  of 
the  variability  in  organic  carbon  change  remained 
unexplained . 

Logarithmic  Transformation 

The  relationship  between  organic  carbon  and  time  is 
frequently  expressed  by  plotting  log  %  of  the  original 
carbon  content3  versus  time  in  order  to  account  for  first 
order  loss.  Log  transformation  did  not  increase  the 
probability  of  a  higher  correlation  (Fig . 26 , 27 , 28  and  29). 

3 

%  of  original  OC  =  (present  %  OC/original  %  OC)100 


Change  in  Concentration  of 
Organic  Carbon  (%) 
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Black  Sampling  Zone 
Continuous  Cropping 


Years  of  Cultivation 


Fig. 25  The  Influence  of  Continuous  Cropping  on  the 

Relationship  Between  Change  in  Concentration  of 
Organic  Carbon  and  Time  in  the  Black  Sampling  Zone 
(r  =  0.38;  p=0.16;  n  =  9 ) . 
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Fig .  26 


Brown  Sampling  Zone 
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7  O 
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.5 
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24  44  65 

Years  of  Cultivation 


Relationship  Between  Log  %  of  Original  Organic 
Carbon  and  Time  in  the  Brown  Sampling  Zone  (r=— 0.51 ; 
p=0 . 06 ;  n=  1  1 ) . 
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Black  Sampling  Zone 


Years  of  Cultivation 


Fig. 27  Relationship  Between  Log  %  of  Original  Organic 

Carbon  and  Time  in  the  Black  Sampling  Zone  (r=-0.05; 
p=0 . 43 ;  n= 1 4 ) . 
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Brown  Sampling  Zone 
Summer  Fallow 


Years  of  Cultivation 


Fig. 28  The  Influence  of  Summerf allowing  on  the  Log  %  of 
Original  Organic  Carbon  and  Time  in  the  Brown 
Sampling  Zone  (r=-0.39;  p=0.17;  n=8). 
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Black  Sampling  Zone 
Continuous  Cropping 


Years  of  Cultivation 


Fig. 29  The  Influence  of  Continuous  Cropping  on  the  Log  %  of 
Original  Organic  Carbon  and  Time  in  the  Black 
Sampling  Zone  (r=0.35;  p=0.18;  n=9). 
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5.3.4  Multiple  Linear  Regression  Models 

Multiple  linear  regressions  were  performed  because  none 
of  the  simple  linear  models  satisfactorily  explained  the 
relationship  between  observed  changes  in  organic  carbon  and 
management  variables.  Three  progressively  more  complex 
models  were  examined  for  each  of  the  dependent  variables 
listed  with  each  set  of  independent  variables. 


First  Model: 


Dependent  Variables: 

present  concentration  of  organic  carbon 
change  in  concentration  of  organic  carbon 
change  in  mass  of  organic  carbon 


Independent  Variables: 

original  concentration  of  organic  carbon 

years  of  cultivation 

average  y ield ( 1 982- 1 967 ) 

crop  rotat ion ( 1 982- 1 967 ) 

organic  residues  returned ( 1 982- 1 967 ) 


Second  Model: 
Dependent  Variables: 

as  before 


Independent  Variables: 

as  before  except  yield,  crop  rotation  and  organic 
residues  were  subdivided  as  follows: 

average  y ield( 1 982- 1 977 )  and  (1977-1967) 
crop  rotat ion ( 1 982- 1 977 )  and  (1977-1967) 
organic  residue  returned ( 1 982- 1 977 )  and  (1977-1967) 
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Third  Model: 
Dependent  Variables: 

as  before 


Independent  Variables: 

as  in  the  second  model  plus 
density  of  Ap 

available  moisture  capacity  in  Ap 
texture  of  parent  material  (hand  textured) 

N  input ( 1982-1977)  and  (1977-1967) 

P  input ( 1982-1977)  and  (1977-1967) 

tillage  intens i ty ( 1 982- 1 977 )  and  (1977-1967) 


On  a  province  level,  the  three  models  consistently 
yielded  the  same  results.  The  change  in  mass  of  organic 
carbon  could  not  be  explained  by  any  of  the  models.  Highest 
prediction  accuracy  (r2=0.6)  for  change  in  concentration  of 
organic  carbon  was  obtained  with  the  second  and  third  models 
which  yielded  identical  results.  Changes  in  concentration  of 
organic  carbon  became  more  negative  with  higher  original 
organic  carbon  content;  but  more  positive  with  continuous 
cropping  and  return  of  organic  residue. 

Y’=  -23.9  -  5.1(original  OC)  +  51.4(crop  rotat ion ( 82-77 )  ) 

-0 . 0005 (yield( 77-67 ) )  +  2.8(organic  residue ( 77-67 )) ; 


n  =  48 . 


'  3 


1  I 
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Original  organic  carbon  explained  26%  of  the  variability; 
21%  were  explained  by  crop  rotation,  9%  by  yield  and  4%  by 
organic  residue.  40%  of  the  variability  remained 
unexplained . 

Yields ( 77-67 )  and  crop  rotat ion ( 82-77 )  explained  43%  of 
the  variability . in  present  concentration  of  organic  carbon: 
Y'=  2.3  -  0 . 0004 (yield ( 77-67 ) )  +  2.6(crop  rotat ion ( 82-77 )) . 

Yield  explained  32%  and  crop  rotation  11%  of  the 
variability.  The  positive  influence  of  continuous  cropping 
on  organic  carbon  content  is  apparent. 

Table  13  summarizes  the  results  obtained  from  the  three 
models  on  a  sampling  zone  level.  It  indicates  variables 
taken  into  the  regression  equations,  the  direction  of  the 
relationship  (positive  slope  or  negative  slope),  and  how 
much  of  the  variability  of  the  dependent  variable  could  be 
explained  by  the  variables  in  the  equation. 

As  at  the  province  level,  prevailing  variables  taken 
into  the  regression  equations  were  management  variables  and 
original  organic  carbon  content.  Within  a  sampling  zone,  a 
high  degree  of  consistency  existed  between  the  three  models 
with  regard  to  the  predictor  variables.  Between  the  sampling 
zones  different  management  variables  seemed  to  be  linked  to 
changes  in  organic  carbon: 

1 .Brown  Sampling  Zone 

Concentration  of  organic  carbon  decreased  with 
increasing  yields  (Table  13;  second  model). 
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Y'  =  -71.4  -  0.009(yield(82~67));  n=11;  1st  model. 

Y’  =  0.98  -  0 . 0003 (yield ( 77-67 )) ;  n  =  8;  2nd  model. 

Change  in  concentration  of  organic  carbon  increased 
with  increasing  yields  (1st  and  2nd  models). 

Y'  =  -78.8  +  0 . 02 (yield( 1 977-67 )) ;  n=8;  2nd  model. 

Return  of  organic  residue  was  positively  linked  with 
mass  changes  in  organic  carbon  and  explained  nearly  60%  of 
the  variability  (2nd  model;  Table  13). 

2 .Dark  Brown  Sampling  Zone 

Present  concentration  of  organic  carbon  increased  with 
increasing  original  concentration  of  organic  carbon  (1st, 2nd 
and  3rd  model).  Prediction  accuracy  more  than  doubled  in  the 
third  model  where  apparent  bulk  density  was  additionally 
taken  into  the  equation. 

Y'  =  0.22  +  0 . 42 ( or iginal  organic  carbon)  -  2.2(Db  cult.); 

n=10;  3rd  model. 

Original  carbon  content  explained  57%,  apparent  bulk 
density  in  the  cultivated  A  horizons  explained  25%  of  the 
variability.  A  strong  inverse  relationship  between  apparent 
bulk  density  and  concentration  of  organic  carbon  is  evident. 
Original  organic  carbon  content  was  moderately  positively 
correlated  with  crop  rotations,  average  annual  yields  and  N 
input  between  1977  and  1967  (Table  14). This  indirectly 
indicates  that  present  concentration  of  organic  carbon  was 
higher  when  more  continuous  cropping  and  higher  N  input  were 
practiced  as  well  as  higher  yields  achieved.  Original  carbon 
content  of  the  soil  may  have  been  related  to  management 
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strategies  employed. 

Only  the  third  model  which  took  soil  properties  apart 
from  concentration  of  organic  carbon  into  account  could 
explain  any  of  the  variability  in  change  in  concentration  of 
organic  carbon.  Apparent  bulk  density  explained  51%  of  the 
variability  in  concentration  change  of  organic  carbon.  The 
change  in  concentration  became  more  negative  when  bulk 
density  increased: 

Yf  =  12.8  -  44.5 (Db  cult.);  n=10. 

Variability  in  the  change  in  mass  of  organic  carbon  was 
explained  consistently  by  crop  rotation,  especially  five 
years  previous  to  sampling  (2nd  and  3rd  model).  The  positive 
direction  of  the  slope  indicates  that  the  change  in  mass 
became  more  positive  the  more  continuous  cropping  was 
practiced . 

Y’  =  -4.6  +  85.5(crop  rotat ion ( 82-77 ))  ;  n=14;  2nd  model. 

Y'  =  -77.8  +  24.5(crop  r ota t i on ( 82- 77 ) ) ;  n-10;  3rd  model. 

3 .Black  Sampling  Zone 

Present  concentration  of  organic  carbon  increased  with 
frequency  of  continuous  cropping  (2nd  model)  and  when  higher 
yields  were  achieved  (3rd  model). 

Changes  in  concentration  of  organic  carbon  became  more 
negative  when  tillage  intensity  increased. 

Y'  =  61.3  -  28.8(tillage  intensity ( 82-77 )) ;  n=9;  3rd  model. 
Tillage  intensity  and  crop  rotation  are  not  independent  from 
each  other  because  tillage  intensity  becomes  higher  the  more 
summer  fallowing  is  practiced.  This  indirectly  indicates  that 
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summer  fallowing  influences  changes  in  concentration  of 
organic  carbon  negatively,  whereas  a  practice  requiring 
fewer  tillage  operations  such  as  continuous  cropping  has  a 
positive  effect. 

The  change  in  mass  of  organic  carbon  became  more 
positive  the  higher  the  original  organic  carbon  content 
(1st, 2nd, 3rd  model);  the  longer  the  cultivation  time  (1st 
and  2nd  model)  and  the  more  continuous  cropping  was 
practiced  (1st  and  2nd  model). 

Y'  =  -153.8  +  1 0 . 6(original  OC)  +  1.2(years  of  cultivation) 

+  98(crop  rotat ion ( 82-67 )) ;  n=13;  1st  model. 

Original  organic  carbon  explained  54%  of  the 
variability,  years  of  cultivation  23%  and  crop  rotation  7%. 
Original  organic  carbon  content  was  positively  correlated 
with  return  of  organic  residue  (Table  14).  This  indirectly 
implies  that  return  of  organic  residue  had  a  positive  effect 
on  changes  in  mass  of  organic  carbon. 

4 .Gray  Luvisolic  Sampling  Zone 

Much  of  the  variability  in  the  change  in  concentration 
of  organic  carbon  seemed  to  be  explained  by  original  organic 
carbon  content  whereby  changes  became  more  negative  the 
higher  the  original  organic  carbon  content  (1st  and  2nd 
model )  . 

Y'  =  -19.3  -  5.1(original  organic  carbon);  n=6;  1st  model. 

Some  of  the  native  profiles  in  this  group  contained 
LFH/Ah  horizons  with  a  high  concentration  of  organic  carbon. 
The  dilution  effect  after  cultivating  became  strongest  in 
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these  cases.  Correlation  of  original  organic  carbon  content 
with  management  variables  as  seen  in  Table  14  and  in  the 
second  model  indicate  that  changes  in  concentrat ion  of 
organic  carbon  became  more  positive  the  more  summer  fallowing 
was  practiced.  Fields  under  summer fallow ,  however,  seemed  to 
have  received  high  inputs  of  mineral  fertilizer  (see  Table  9 
and  10). 

The  second  model  indicated  on  the  other  hand,  that  the 
change  in  mass  became  more  positive  when  more  organic 
residue  was  returned.  High  values  for  organic  residue  coded 
in  the  Gray  Luvisolic  area  were  for  input  of  manure  and 
green  manure  in  combination  with  lay  farming.  A  positive 
change  of  91%  in  mass  of  organic  carbon  in  the  A  horizon  was 
found  for  one  site  under  lay  rotation  (Appendix  II. 5). 
Further  interpretation  of  these  data,  however,  is  restricted 
due  to  the  low  number  of  cases  in  the  model. 

6 .Peace  River  Sampling  Zone 

Close  to  60%  of  the  variability  in  the  change  in 
concentration  of  organic  carbon  was  accounted  for  by 
original  organic  carbon  content.  Changes  became  more 
negative  for  sites  where  carbon  content  was  originally 
higher.  This  was  the  case  for  the  Gray  Luvisolic  soils  in 
the  Peace  River  area.  Correlation  of  original  organic  carbon 
content  with  management  variables  (Table  14)  implies  that 
losses  in  concentration  of  organic  carbon  became  higher  when 
more  organic  residue  was  returned.  Highest  return  of  organic 
residue  was,  however,  to  one  of  the  Gray  Luvisolic  sites; 


data  for  the  other  site  were  missing.  Considering  the 
variety  of  soils  in  this  area  the  number  of  cases  is  low  and 
restricts  interpretation  and  conclusions. 

The  third  model  indicated  that  P  input  in  the  five 
years  prior  to  sampling  had  a  positive  effect  on  changes  in 
concentration  of  organic  carbon.  Highest  P  input  was  to 
fields  under  continuous  cropping  (Table  10).  Greater  losses 
seemed  to  have  occurred  on  soils  with  finer  textures  (3rd 
model) (for  texture  codes  see  Appendix  I,  codebook  farminfo). 
Y'  =  -15  +  4.2(P(82-77) )  -  1 5 . 2 ( org . residue ( 77-67 ) ) 

-  6.0(texture  of  parent  material);  n=7;  3rd  model. 

Additions  of  P  explained  74%  of  the  variability, 
organic  residue  an  additional  20%,  and  texture  of  the  parent 
material  an  additional  5%. 

Changes  in  mass  of  organic  carbon  became  more  negative 
when  higher  yields  were  achieved  (1st, 2nd, 3rd  model). 

Y'  =  147.8  -  0.05(yield(82-77) ) ;  n=7;  2nd  model. 

They  became  more  negative  with  higher  input  of  N  and 
organic  residue  (3rd  model),  but  more  positive  when  more  P 
was  applied  (3rd  model). 

Y'  =  174.6  -  0.06(yield(82-77) )  -  0 . 45 (N ( 82-77 ) ) 

+  1 . 26 ( P ( 77-67 ) )  -  1 . 3 ( org . residue ( 82-77 )) ;  n=7; 

3rd  model. 

Yield  explained  85%  of  the  variability,  14%  was 
explained  by  N,  0.9%  by  P,  and  0.1%  by  organic  residue. 

Interpretation  is  limited  by  the  low  number  of  cases 
together  with  soils  of  different  Great  Groups  in  this  area. 


5.4  Discussion 


Management  has  changed  since  the  late  1960s.  This  was 
mainly  a  reduction  in  summer  fallowing  and  an  increase  in 
continuous  cropping  with  a  concomitant  increase  in  return  of 
organic  residue  and  in  mineral  N  and  P  input.  This  seemed  to 
have  been  associated  with  an  increase  in  yield.  This  is 
consistent  with  general  trends  in  the  prairie  provinces  as 
documented  by  Ridley  and  Hedlin  (1980).  The  authors  found  a 
positive  correlation  between  increased  fertilizer  input  and 
yields.  Further  reasons  were  biological  and  technological 
progress  which  is  difficult  to  assess.  The  historical  trend 
reported  by  Ridley  and  Hedlin(1980)  of  decreasing  yields  due 
to  drought  and  erosion  in  the  dirty  thirties  with  a 
succeeding  increase  in  yields  is  not  evident  from  the  data 
of  this  study.  A  more  thorough  survey  necessary  for  such  a 
purpose  was  not  within  the  scope  of  this  thesis. 

Average  annual  cereal  yields  in  the  15  years  previous 
to  sampling  were  1 . 6 ( i  f  the  one  farm  with  continuous 
cropping  is  eliminated),  1.7,  2.3,  1.8  and  2.1  Mg  ha-1  in 
the  Brown,  Dark  Brown,  Black,  Gray  Luvisolic  and  Peace  River 
zones  (Table  11).  They  were  1.6,  2.3,  1.7  and  2.2  Mg  ha'1 
for  the  Dark  Brown,  Black,  Gray  Luvisolic  and  Peace  River 
zones  between  1977  and  1967  (Table  11).  Average  cereal 
yields  in  Alberta  based  on  1966-1974  Alberta  Hail  and  Crop 
Insurance  Corporation  records  were  1.7,  2.3,  2.4,  1.6  and 
1.5  Mg  ha"1  in  the  Brown,  Dark  Brown,  Black,  Gray  Luvisolic 
and  Peace  River  zones.  These  data  were  calculated  from 
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Tables  D1  to  D15  of  Peters  and  Pettapiece ( 1 98 1 ) .  Data  of  the 
present  study  seem  to  underestimate  yields  in  the  Dark  Brown 
zone  and  overestimate  yields  in  the  Gray  Luvisolic  and  Peace 
River  zone.  A  comparatively  low  number  of  cases  was 
available  in  this  study,  but  data  do  not  diverge  markedly 
from  provincial  averages  obtained  from  Peters  and 
Pet tapiece ( 1 98 1 )  using  a  much  larger,  data  base.  These  data 
suggest  the  sampling  of  farm  units  was  not  biased  toward  any 
group  of  farmers. 

The  change  to  continuous  cropping  affected  more  the 
central  and  northern  regions  of  the  province  as  indicated  by 
the  present  percentage  of  summer fallowing  in  the  various 
sampling  zones.  Similar  trends  were  reported  by  Anderson  & 
Associates  Ltd. (1981).  Factors  explaining  summerfallow 
frequency  in  their  study  were  mainly  debt-equity  ratios  and 
land-factor  cost  ratios.  Pre-season  moisture,  weed  control 
and  energy  conservation  played  a  subordinate  role.  The 
authors  indicated  that  summer  fallowing  could  be  a  very 
rational  practice  from  a  private  point-of -view  on  coarser 
drier  soil  in  southern  regions  of  the  province.  They  further 
concluded  that  economic  benefits  from  summer fal low ing  were 
short-run  and  accrued  directly  to  farmers  whereas  many  of 
the  costs  were  long-run  and  accrued  to  society.  Such 
soc io-economic  costs  are  salinity  problems,  erosion  and, 
linked  to  it,  organic  matter  loss. 

Statistical  analysis  at  the  sampling  zone  level  yielded 
a  higher  resolution  than  was  possible  at  the  province  level, 


indicating  regional  aspects  are  important. 

None  of  the  simple  linear  models  fully  explained  the 
variability  in  present  concentration  of  organic  carbon  and 
in  changes  in  concentration  and  mass  of  organic  carbon  with 
time.  This  was  not  unexpected  as  the  present  models  on 
organic  carbon  dynamics  suggest  a  continuous  decline  in 
organic  matter  with  time  according  to  a  first  order 
reaction . 


A  curvilinear  trend  for  declining  soil  organic  carbon 
with  time  which  is  in  agreement  with  existing  models 
appeared  to  be  present  only  for  the  change  in  concentration 
of  organic  carbon  in  the  Brown  sampling  zone.  For  the 
remaining  zones,  as  well  as  for  the  whole  province,  there 
was  no  evidence  of  a  consistent  continuing  decline  in 
organic  carbon  with  time,  but  some  evidence  of  a  positive 
trend.  This  suggests  that  although  organic  carbon  has 
decreased  since  cultivation  started,  as  shown  in  the 
previous  chapter,  there  might  be  an  increasing  trend  at 
present.  This  was  supported  by  results  from  loglinear 
tranf ormat ions ,  and  is  consistent  with  changes  in  cropping 
systems  during  this  century. 

Multiple  linear  regression  models  indicated  that 
positive  changes  in  concentration  and  mass  of  organic  carbon 
were  related  to  higher  yields  and  return  of  organic  residue 
in  the  Brown  sampling  zone;  they  were  related  to  continuous 
cropping  in  the  Dark  Brown  sampling  zone.  Mass  changes  in 


the  Black  sampling  zone  were  positively  related  to  original 
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concentration  of  organic  carbon,  years  of  cultivation  and 
continuous  cropping.  In  the  Gray  Luvisolic  and  Peace  River 
regions,  however,  changes  in  concentration  and  mass  were 
negatively  influenced  by  original  concentration  of  organic 
carbon.  This  shows  the  effect  of  the  rapid  loss  from  LFH 
material  on  the  one  hand  and  slower  loss  of  Ah  material  on 
the  other.  Such  observations  are  consistent  with  the 
hypothesis  that  association  with  mineral  particles  in  the  Ah 
horizon  provides  physical  protection  to  organic  carbon. 

Whereas  continuous  cropping,  higher  yields  and 
increased  return  of  organic  residue  had  positive  effects  on 
changes  in  organic  carbon  in  the  Chernozemic  zones,  it 
seemed  to  have  the  opposite  effect  in  the  Peace  River  and 
Gray  Luvisolic  zones.  Although  interpretation  is  restricted 
by  a  low  number  of  cases  in  the  Luvisolic  zone  and  a  range 
of  Great  Groups  in  the  Peace  River  zone,  this  indicates  that 
the  original  nature  of  the  soil  determines  its  reaction  to 
management  practices.  The  rapid  loss  from  LFH  material  seems 
to  have  an  overriding  effect  on  ameliorative  practices  like 
return  of  organic  residue  in  the  Peace  River  zone.  The 
residues  were  added  to  the  site  most  susceptible  to  loss  and 
dilution.  The  loss  of  organic  carbon  with  more  residues  in 
the  Peace  River  area  was  therefore  only  a  correlation  and 
not  an  effect.  Linear  regression  models  do  not  describe 
these  complex  relationships. 

On  a  provincial  level,  the  change  in  concentration  of 
organic  carbon  was  negatively  influenced  by  original  organic 
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carbon  content  and  yield;  this  suggests  the  influence  of  the 
Peace  River  area.  They  were  positively  affected  by 
continuous  cropping  and  organic  residue  indicating  an 
influence  of  the  Chernozemic  zones.  The  negative  effect  of 
yield  is  likely  related  to  the  fact  that  yield  was  reported 
and  used  for  the  year  obtained  and  not  as  annual  yield. 
Consequently,  fallow  systems  may  have  highest  yield  in  the 
non-fallow  year  but  lowest  annual  yield,  which  then  shows  up 
as  a  negative  yield  effect.  Therefore,  such  methods  of 
increasing  yields  are  destructive  of  organic  matter. 

Table  14  suggests  a  relationship  between  original 
organic  carbon  content  and  management  practices  employed. 
This  relationship  might  reflect  the  capability  of  the  soil 
for  crop  production.  Similar  trends  were  reported  by 
Reinl ( 1 98 1 )  .'  Original  concentration  of  organic  carbon 
reflects  a  specific  combination  of  soil  forming  factors  such 
as  climate  or  parent  material  which  are  likely  to  influence 
management  decisions. 

Statistical  analysis  in  this  section  dealt  only  with 
data  from  the  last  15  years  prior  to  sampling.  Too  many  data 
from  1967  to  1907  were  missing  to  allow  statistical 
analysis.  The  data  base  indicates,  however,  that  changes  in 
management  occurred  mainly  during  these  last  15  years 
previous  to  sampling.  The  main  changes  are  a  substitution  of 
summerfallow  rotations  by  continuous  cropping  with  a 
concomitant  increase  in  the  amount  of  input  of  mineral 
fertilizer  and  straw  return.  This  might  be  linked  to  the 
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observation  that  no  consistent  continuing  decline  in  organic 
carbon  over  time  was  evident  and  to  the  hypothesis  of  a 
trend  toward  increasing  organic  matter  following  an  initial 
decline.  The  fact  that  higher  losses  occurred  under 
summerfallow  than  under  continuous  cropping  is  consistent 
with  the  above  hypothesis.  Rotation  studies  have  shown  that 
continuous  cropping  is  less  destructive  to  soil  organic 
matter  than  is  summer  fallowing .  This  is  due  mainly  to  the 
increased  input  of  belowground  residues.  Steady-state  levels 
of  organic  matter  are  determined  by  the  ratio  of  addition  to 
removal  rates.  An  increase  in  addition  rate  will  cause  a 
proportional  increase  in  the  steady-state  quantity  if 
removal  rate  remains  constant.  A  change  to  a  management 
system  with  a  higher  ratio  of  addition  rate  to  removal  rate 
will  therefore  cause  organic  carbon  to  increase  toward  a 
higher  steady-state  level.  Existing  mathematical  models  of 
organic  matter  dynamics  treat  management  as  constant  over 
time.  The  present  study  has  shown  that  inputs  of  organic 
residue  have  increased  since  cultivation  started  and  are  not 
constant  with  time.  Therefore,  the  addition  term  (R)  in 
these  models  must  be  modified  to  describe  the  present 
situation . 


5.5  Conclusions 

Management  practices  changed 
Since  then  continuous  cropping 
straw  return  and  input  of  N 


since  the  late  1960s. 
with  concomitantly  higher 
and  P  has  replaced 


summer fallowing  especially  in  the  central  and  northern 
regions  of  the  province.  Increases  in  yield  took  place  in 
the  same  period  especially  under  continuous  cropping.  Losses 
in  organic  carbon  tended  to  be  lower  under  continuous 
cropping  than  under  summer fallow . 

No  simple  linear  model  could  explain  the  variability  in 
organic  carbon  changes  with  time.  A  higher  level  of 
resolution  was  reached  at  the  regional  than  at  the 
provincial  level. 

There  is  no  evidence  of  a  consistent  continuing  decline 
of  organic  carbon  with  time.  This  holds  for  the  whole 
province  as  well  as  for  the  sampling  zones.  Only  in  the 

f 

Brown  sampling  zone,  was  a  negative  curvilinear  trend 
evident  for  changes  in  concentration  of  organic  carbon  with 
time.  In  the  remaining  sampling  zones  no  such  trend  was 
obvious.  This  was  supported  by  results  from  loglinear 
transformations . 

Multiple  linear  regression  models  indicated  that  the 
change  in  organic  carbon  was  positively  related  to 
continuous  cropping,  yields  and  original  organic  carbon 
content  in  the  Brown,  Dark  Brown  and  Black  sampling  zones. 
Interpretation  was  limited  for  the  Gray  Luvisolic  and  Peace 
River  area  due  to  low  number  of  cases.  The  linear  treatments 
indicated,  however,  that  the  effect  of  continuous  cropping, 
organic  residue  and  yields  in  the  Chernozemic  zones  was 
opposite  to  that  in  the  Luvisolic  and  Peace  River  zones. 
This  seemed  to  be  related  to  the  original  nature  of  the  soil 
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in  these  areas.  Loss  of  organic  carbon  with  ameliorative 
practices  was  therefore  a  correlation,  not  an  effect. 

Changes  in  management  in  the  last  15  years  previous  to 
sampling  seem  to  have  reversed  the  direction  of  organic 
carbon  changes  toward  a  higher  steady-state  level. 


6.  Overall  Summary  and  Conclusions 

Significant  changes  in  organic  carbon  have  resulted  from 
cultivation  in  Alberta.  On  the  sampling  zone  level,  changes 
were  negative.  Concentration  of  organic  carbon  in  cultivated 
A  horizons  decreased  by  52%,  42%,  40%,  13%,  75%  and  64%  in 
the  Brown,  Dark  Brown,  Black,  Dark  Gray,  Gray  Luvisolic  and 
Peace  River  sampling  zones.  Mass  in  the  A  horizon  decreased 
by  31%,  28%,  17%,  29%,  27%  and  15%  for  the  same  sampling 
zones.  Mass  in  the  profile  decreased  by  36%,  25%,  21%,  22%, 
39%  and  28%  in  these  sampling  zones.  At  the  individual  farm 
level,  however,  some  positive  changes  in  mass  of  organic 
carbon  occurred:  21%,  13%,  33%,  43%  and  42%  of  the  sites  in 
the  Brown,  Dark  Brown,  Black,  Gray  Luvisolic  and  Peace  River 
sampling  zones  had  a  gain  in  mass  of  organic  carbon  in  the  A 
horizon.  A  mass  gain  in  the  profile  occurred  on  14%,  13%, 
14%  and  33%  of  the  sites  in  the  Brown,  Black,  Gray  Luvisolic 
and  Peace  River  area.  These  farms  were  viable  units.  One 
must  conclude , therefore ,  that  it  is  possible  to  increase 
organic  matter  and  maintain  a  viable  farm  operation.  The 
financial  aspect  of  this  merits  further  detailed 
i nvest igat ion . 

The  estimate  of  McGill  et  al(l981)  of  a  35%-50%  organic 
matter  loss  on  the  Canadian  Prairies  lies  within  the  range 
of  concentration  and  mass  losses  but  tends  to  be  higher  than 
those  found  here.  Data  reported  in  the  literature  on  organic 
matter  loss  correspond  with  14%  of  the  farms  in  the  Brown 
and  17%  of  those  in  the  Dark  Brown  sampling  zones  in  the 
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present  study  under  similar  crop  rotation  and  cultivation 
time.  Most  literature  values  are  derived  from  only  one  to 
three  sites,  however.  They  can  not  appropriately  represent 
the  average  for  any  soil  zone.  A  larger  number  is  needed 
than  has  been  used  in  the  past. 

Loss  in  concentration  tended  to  be  higher  than  losses 
in  mass.  This  can  be  attributed  to  dilution  resulting  from 
subsurface  horizon  material  being  mixed  into  the  A  horizon. 
The  thickness  of  A  horizons  increased  in  cultivated 
profiles,  though  not  significantly;  B  horizon  thickness 
tended  to  decrease.  Linked  to  this  was  an  increase  in 
apparent  bulk  density.  An  inverse  relationship  was  found 
between  changes  in  apparent  bulk  density  and  changes  in 
concentration  of  organic  carbon. 

Moisture  retention  at  -33  and  -1500kPa  in  A  horizons 
decreased  due  to  cultivation.  The  effect  was  greater  at 
lower  potentials  than  at  higher  potentials  with  the  result 
that  available  moisture  capacity  increased  due  to 
cultivation . 

No  reduction  in  profile  depths  occurred  upon 
cultivation  suggesting  no  major  erosion  effects  on  organic 
matter  losses.  Variability  in  profile  thicknesses  causes 
such  a  measure  of  erosion  to  have  a  low  resolution,  however, 
and  a  definitive  answer  to  the  relation  between  erosion  and 
decomposition  awaits  further  study. 

Simple  linear  regressions  and  log  transformations 
indicated  that  there  was  no  consistent  decline  in  organic 
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carbon  with  time  over  the  whole  province  nor  in  individual 
sampling  zones.  Only  the  change  in  concentration  of  organic 
carbon  with  time  in  the  Brown  sampling  zone  showed  the 
negative  curvilinear  trend  embodied  in  existing  concepts. 
There  is  a  possibility  that  losses  have  ceased.  This  seemed 
to  be  linked  to  changes  in  management  practices  during  the 
last  5  to  15  years  previous  to  sampling  in  1982.  During  that 
time  continuous  cropping  with  more  straw  return  and 
application  of  mineral  fertilizer  became  more  frequent. 
Yields  increased  simultaneously,  mainly  under  continuous 
cropping . 

Multiple  linear  regression  analyses  suggested  that 
present  higher  concentrations  of  organic  carbon  and  changes 
in  concentration  and  mass  of  organic  carbon  were  positively 
influenced  by  continuous  cropping  and  yields  in  the  Brown, 
Dark  Brown  and  Black  sampling  zones.  The  lowest  loss  in 
organic  carbon  in  the  Black  sampling  zone  therefore  seems  to 
be  related  to  the  highest  frequency  of  continuous  cropping 
and  highest  average  yields  in  the  last  15  years  previous  to 
sampling.  The  Brown  sampling  zone  on  the  other  hand  had  the 
highest  proportion  of  summerfallow  and  highest  proportional 
loss  of  organic  carbon  of  the  Chernozemic  zones. 
Interpretation  of  data  for  the  Gray  Luvisolic  and  Peace 
River  area  was  limited  due  to  small  sample  size.  The  models 
indicated,  however,  that  the  soils  in  these  regions  reacted 
opposite  to  the  soils  in  the  Chernozemic  areas  when  under 
the  same  management.  This  might  be  related  to  the  original 
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nature  of  the  soil.  It  shows  the  effect  of  the  rapid  loss 
from  LFH  material  as  opposed  to  physical  protection  of 
organic  carbon  when  associated  with  mineral  particles  in  Ah 
horizons.  This  rapid  loss  had  an  overriding  effect  on 
ameliorative  practices  like  returning  of  organic  residues. 

Results  for  the  Brown,  Dark  Brown  and  Black  zones 
suggest  that  change  in  management  to  continuous  cropping 
together  with  higher  input  of  mineral  fertilizer  has 
reversed  the  direction  of  changes  in  organic  carbon.  This  is 
consistent  with  the  literature  from  rotation  studies.  It  has 
never  before  been  shown  at  the  farm  level.  Losses  of  organic 
carbon  under  continuous  cropping  were  lower  than  under 
summerfallow  in  the  Brown,  Dark  Brown  and  Black  sampling 
zones . 

Some  of  the  above  trends  were  not  statistically 
significant,  however.  This  could  be  due  to  insufficient 
sample  size  per  sampling  zone  even  though  this  study  is  the 
most  thorough  yet  reported  in  the  literature.  Data  obtained 
are  the  first  data  base  derived  from  farm  sites  on  a  broad 
geographic,  climatic  ,  soil  and  management  base.  Existing 
data  which  are  the  basis  for  present  theoretical  models  were 
derived  from  few  sample  sites,  research  plots  or  lab 
experiments.  Data  of  this  study  may  be  used  in  testing 
existing  theoretical  models.  These  models  will  have  to  be 
reviewed  to  provide  a  linkage  for  the  above  findings. 
Mathematical  models  have  up  to  now  treated  management  as  a 
constant  and  organic  carbon  decreased  continuously  in  a 


124 


series  of  first  order  reactions  to  lower  levels  than  found 
here.  Results  of  this  thesis  have  shown  that  changes  in 
management  do  occur  and  can  reverse  changes  in  organic 
carbon.  Mathematical  models  should  therefore  be  modified  to 
accommodate  such  changes  in  management  and  their  effects. 
Modifications  with  regard  to  changing  inputs  of  organic 
residues  with  time  would  make  the  present  models  more 
consistent  with  the  results  of  the  present  study. 


' 
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7.  APPENDIX  I 


Sample  Ques 
Codebook  fo 
Codebook  fo 
Codebook  fo 


ionna ire 
File  Farminfo 
File  Cropinfo 

File  Analysis:  Labanalysis  Data 


/ 
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QUESTIONNAIRE 


Name : 

Address : 

Legal  Location: 

Soil  Zone: 

Soil  Classification: 

Total  time  cultivated  years: 


Cropping  systems  since  beginning  of  cultivation : 
Year  Cropping  system 


Application  of  chemical  fertilizer : 

Year  Amount  of  N  P  K  S  Lime 


Organic  residue  returned : 


Year 

Straw  removed 
Stubble  Mulch 
All  Straw  Returned 
Manure 

Green  Manure 


13  0  .  nuz,:di 
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Tillage  operations  during : 


Cropping  systems: 


Fallow : 


Weed  Control: 


chemical 

mechanical 

cultural 


Yields : 


Year 

Wheat 

Barley 

Oats 

Canola 


How  would  you  rate  erosion  effects  since  cultivation  started? 

severe? .  light? . 

moderate? .  negligible?.... . 

Dominant  type  of  erosion:  wind .  water . 

Did  you  notice  any  changes  in  organic  matter  content? 

No .  If  yes,  which?  . . . 

Did  you  notice  any  changes  in: 

Tilth:  No  .  If  yes,  which?  . . . 

Water  Storage:  No  .  If  yes,  which?  . 

Nutrient  Supply:  No  .  If  yes,  which?  . 

What  are  your  future  management  plans?  . . . 
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Codebook  for  File  Farminfo 
Col  1-  2:  Farm- ID  01  -  80 


Col  4:  Land  use  1=cultivated 

2=nat i ve 


Col  6-14:  Legal  location  1 =NW  2=NE  3=SE  4=SW  5=centre 

Sect ion=0 1 -36 
Township=01 - 120 
Range=0 1 -30 
Meridian  4=W4  5=W5 


Col  16:  Soil  order  1=Brunisolic 

2=Chernozemi c 
3=Gleysol ic 
4=Luvi sol ic 
5=0rgani c 
6= Podzol ic 
7=Regosol ic 
8=Solonetzi c 
9=undef i ned 


Col  18-19:  Great  group  1=Melanic  Brunisol 

2=Eutric  Brunisol 
3=Sombric  Brunisol 
4=Dystric  Brunisol 
5= Brown 
6=Dark  Brown 
7  =  B 1 ack 
8=Dark  Gray 
9=Humic  Gleysol 
1 0=Gleysol 
1 1 =Luvi c  Gleysol 
12=Gray  Luvisol 
1 3=F ibr i sol 
1 4=Mesi sol 
1 5=Humi sol 

16=Humo  Ferric  Podzol 

1 7=Regosol 

18=Humic  Regosol 

1 9=Solonetz 

20=So Iodized  Solonetz 

2 1 =Solod 

99=undef i ned 


Col  21-24:  Subgroup 


1=alkal ine 
2=black 


. 
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3=brown 
4=bruni sol ic 
5=ca lcareous 
6=cumul ic 
7=darK 

8=darK  brown 
9=darK  gray 
1 0=eluvi ated 
1 1  =  f err i c 
12=fera 
13=gleyed 
1 4=gray 
1 5=humi c 
1 6=or  thi c 
17=podzol ic 
1 8=rego 
1 9=solonetzi c 
20= typic 
99=undef i ned 


Col  26-28:  Series  999=undef i ned 

Col  30-31 :  Ecoregion  99=undefined 


Col  33-34:  Physiographic  section  99=undefined 


Col  36:  Agroclimate 


1 =Rai nf a  1 1 
frost  -  free 
2=Rai nf al 1 
vegetation 
frost  -  free 

3  =  Rai nf a  1 1 
frost  -  free 

4  =  Rai nf a  1 1 
days  frost- 
5=Rainfal 1 
days  frost- 
6  =  Rai nf a  1 1 
frost-free 
7=Mountai n 
precambr i an 


adequate; >90  days 

limiting  50%  of 
time;-  >90  days 

limiting;  >100  days 

adequate;  75-90 
free 

adequate;  60-75 
free 

adequate;  <60  days 

complex ; 
shield 


Col  38:  Soil  Capability  1=land  has  no  significant 

limitations  in  use  for  crops 
adapted  to  area 
2,3=moderate  to  moderately 
severe  limitations  that 
restrict  range  in  crops  or 
require  special  conservation 
practices;  land  is  good  to 


■ 
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fair  for  agricultural 
production 

4=severe  limitations,  range  of 
crops  restricted,  special 
conservation  practices; 
marginal  for  sustained  cereal 
crop  production;  low 
productivity  for  wide  range  of 
crops  or  good  productivity  for 
restricted  range 
5,6=severe  limitations, 
restricted  capability  to 
production  of  perennial  forage 
crops  or  native  pasture 
7=no  capability  for  arable 
agricultural  or  permanent 
pasture 
9=undef i ned 


Col  40-41:  Landform  1=apron 

2=de 1 ta 
3=  fan 

4= 1 ake  plain 
5= hummocky 
6= i nc 1 i ned 
7= level 
8=r idged 
9=rol 1 i ng 
1 0  =  terraced 
1 1 =undu 1  at i ng 
99=undef i ned 


Col  43:  Parent  Material 


1 =g 1 aci of  1 uvi a  1 
2=g 1 aciol acustr i ne 
3=aeol i an 
4=f luvi al 
5= 1 acustr i ne 
6=col luvi al 
7=organic 
8=  t  i  1  1 
9=undef ined 


Col  45:  Texture  of  parent  material 

1=coarse  sand 
2=medium  sandy 
3=fine  sandy 
4=sandy  loam 
5=s i 1 ty  loam 
6=clayey  loam 
7  =  si  1  ty 
8=clayey 


9=organic 
0=undef i ned 


Col  47-53:  Color  of  A  horizon  Hue:  1=1 OYR  2=2. 5YR 

3  =  7. 5YR  4  =  5Y  5  =  5G 
Value:  as  determined 
Chroma:  as  determined 


Col  55-56 :  Years  of  cultivation  as  determined 

Col  58:  Erosion  type  1=wind 

2=water 

3=both 

0=no  erosion  Col  59-60 

Erosion  extent 

col  59 :  overa 1 1 

col  60 :  hill  top 

1 =severe 

2=moderate 

3= 1 ight 

4=neg 1 igible 

0=no  erosion 


Col  62:  Noticed  changes  in  organic  matter 

1 =strong  loss 
2=med i um  loss 
3=light  loss 
4=gai n 

5  =  no  changes 
0=nat i ve 


Col  64:  Noticed  changes  in  tilth 

1  =no 

2=yes , i mproved 
3= yes .deteriorated 


Col  66:  Noticed  changes  in  water  storage 

1  =no 

2=yes , improved 
3= yes .deter i or a ted 


Col  68:  Noticed  changes  in  nutrient  supply 

1  =no 

2=yes , i mproved 
3= yes .deter i or a ted 


Col  70-75:  Future  Management  01=more  fertilizer 


V 
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CodebooK  for  File  Cropinfo 

File  cropinfo  contains  5  records  per  farm; the  5  lines 
represent  different  years: 


1 

2 

3 

4 

5 


1982-1978 

1977-1967 

1966-1947 

1946-1927 

1926-1907 


Col  1-  2:  farm  ID 


Col  4:  crop  rotation:  1 =crop-summer f a  1  low; 

2 =crop- crop -summer fa  1  low; 3=cont i nuous 
croppi ng ; 4=cont i nuous  pasture;  5=lay 
f armi ng ; 6=nat i ve ; 9=mi ssi ng 


Col  6:  main  crop  in  the  system:  1=cereal;  2=forage; 

3=oi 1  crop;  4=legume;  6=native;  9=missing 


Col  8-11:  average  yield  of  main  crop  for  specified 
yearsTkg/ha ) 


Col  13:  second  crop  in  system 


Col  15-18:  average  yield  of  second  crop(kg/ha) 


Col  20:  organic  residue  returned;  1=straw  baled;  2=straw 
burnt  3=stubble  mulch;  4=manure;  5=green  manure; 
6=a 1 1  straw  returned;  7=native;  9=missing 


Col  22-25:  N  fertilizer  (kg/ha) 


Col  27-30:  P  (kg/ha) 


Col  32-35:  K  (kg/ha) 


Col  37-40:  S  (kg/ha) 


Col 


42:  Tillage  intensity:  1=Zero;  2=minimum; 


,  I 
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3=convent iona 1 ;  4=excessive;  6=native; 


Col 


44 


Weed  control:  1 =chemi ca 1 ;  2=mechanical ; 
3=cultural;  6=native;  9=missing 


9=mi ssi ng 


/ 


. 
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Codebook  for  File  Analysis:  Labanalysis  Data 


Col 

1-  2: 

Farm- ID 

0  -  80 

Col 

4: 

Horizon 

1  = Ap/ Ah 

2  =  Ae 

3  =  AB 

4  =  B 

5  =  BC 

6  =  C 

Col 

6-  7: 

Repl icate 

01 -05=cul tivated 
06- 1 0=nat i ve 

Col 

9-15: 

Weight (g) 

Col 

17-20: 

Depth(cm) 

Col 

22-25: 

Compact  ion (cm) 

Col 

27-29: 

pH 

Col 

31-33: 

Inorganic  carbon(%) 

Col 

36-39: 

Total  carbon(%) 

Col 

41-44: 

Organic  carbon(%) 

Col 

46: 

Slope  posi tion 

1  =crest 

2=upper 
3=mid 
4= lower 
5=  toe 

6=depression 
9=undef i ned 
0=no  slope 


Col 

48-51 : 

Moisture 

retention ( % ) 

at 

1/3  bar 

Col 

53-56: 

Moisture 

retention ( % ) 

at 

1 5  bar 

Col 

58-61  : 

%  available  moisture 

. 


8.  APPENDIX  II 


Summary  Results  of  Cultivated  and  Native  Sites 
by  Farm  in  the  Sampling  Zones. 

Each  row  of  data  represents  a  separate  farm. 

*  indicates  differences  at  the  5%  significance  level  between 
cultivated  and  native  sites. 

**  indicates  differences  at  the  1%  significance  level 
between  cultivated  and  native  sites. 
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Appendix  II. 1  BROWN  SAMPLING  ZONE 
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Appendix  III.4:  Simple  Linear  Regression  Analysis  for  the  PROVINCE:  Influence  of  Management  Pract 1 ces (  1 982-  1 967  )  on  the 

Relationship  Between  Present  Concentration  of  Organic  Carbon(l),  Change  In  Concentration  of  Organic 
Carbon(2)  and  Change  In  Mass  of  Organic  Carbon(3)  with  Cultivation  Time.  . 
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Appendix  III. 5:  Simple  Linear  Regression  Analysis  for  the  BROWN  SAMPLING  ZONE:  Influence  of  Management 

Pract1ces( 1982-1967)  on  the  Relationship  Between  Present  Concentration  of  Organic  Carbon(l),  Change  In 
Concentration  of  Organic  Carbon(2)  and  Change  In  Mass  of  Organic  Carbon(3)  with  Cultivation  Time. 
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Append i x  1 1 1 . 6  :  Simple  Linear  Regression  Analysis  for  the  DARK  BROWN  SAMPLING  ZONE:  Influence  of  Management 

Pract 1ces(  1982- 1967)  on  the  Relationship  Between  Present  Concentration  of  Organic  Carbon(l),  Change  In 
Concentration  of  Organic  Carbon(2)  and  Change  in  Mass  of  Organic  Carbon(3)  with  Cultivation  Time. 
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Appendix  III  .7:  Simple  Linear  Regression  Analysis  for  the  BLACK  SAMPLING  ZONE:  Influence  of  Management 

Pract1ces( 1982-1967)  on  the  Relationship  Between  Present  Concentration  of  Organic  Carbon(l),  Change  in 
Concentration  of  Organic  Carbon(2)  and  Change  in  Mass  of  Organic  Carbon(3)  with  Cultivation  Time. 
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Appendix  III .8:  Simple  Linear  Regression  Analysis  for  the  GRAY  LUVISOLIC  SAMPLING  ZONE:  Influence  of  Management 

Pract 1ces( 1982- 1967)  on  the  Relationship  Between  Present  Concentration  of  Organic  Carbon(l),  Change  In 
Concentration  of  Organic  Carbon(2)  and  Change  In  Mass  of  Organic  Carbon(3)  with  Cultivation  Time. 
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DV :  dependant  variable  1=present  concentration  of  organic  carbon;  2=change  In  concentration  of  organic  carbon 

3=change  In  mass  of  organic  carbon 
r:  correlation  coefficient;  Sign.:  significance; 

Y-Int.:  Intercept  with  Y-axis;  N:  number  of  cases 


